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FOREWORD 


This  report  was  prepared  by  Cornell  Aeronautical  Laboratory,  Inc.  , 
Buffalo,  New  York  under  USAF  Contract  No.  AF  33(6l5)-67-C- 1 144,  Project 
No.  651-  E.  The  work  was  administered  under  the  direction  of  the  Air  Force 
Aero  Propulsion  Laboratory,  Research  and  Technology  Division,  Mr.  Harold 
F.  Chambers,  Jr.,  Contract  Monitor. 

This  is  the  project  final  technical  report  that  describes  and  summarizes 
studies  conducted  during  the  period  from  April  1967  through  September  1967 
on  analog  and  digital  data  reduction  techniques  for  heat  transfer  measurements 
made  with  thin-film  resistance  thermometers.  A  draft  of  this  report,  iden¬ 
tified  by  the  contractor  as  Report  No.  AE-2385-Y-2,  was  submitted  on 
1  October  1967. 

This  effort  was  performed  jointly  by  the  Hypersonic  Facilities  Depart¬ 
ment,  the  Computer  Services  Department,  and  the  Avionics  Department  of 
Cornell  Aeronautical  Laboratory  (CAL).  Acknowledgement  is  made  of  the 
major  technical  contributions  of  the  following  individuals:  V.  V.  Abrahamian, 
analog  computer  (COMCOR)  studies,  W.  D.  Fryer,  analog  computer  studies 
and  digital  program  development,  D.  H.  Bock,  P.  F.  Heggs,  and  S.  A.  R.  Ayyar 
analog  circuit  development. 

Publication  of  this  report  does  not  constitute  Air  Force  approval  of  the 
report's  findings  or  conclusions.  It  is  published  only  for  the  exchange  and 
stimulation  of  ideas. 
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ABSTRACT 


An  account  is  presented  o£  experimental  studies  of  techniques  of  data 
reduction  required  to  convert  measurements  obtained  with  thin-film  resistance 
thermometers  to  incident  heating  rates.  Emphasis  is  on  development  of 
practical  data  reduction  techniques  accommodating  nonspecific  gage  surface 
temperature  histories  as  well  as  temperature  excursions  sufficiently  large 
to  necessitate  correction  for  variable  thermal  properties  of  the  gage. 

Analog  computer  solutions  were  obtained  for  a  range  of  surface  tempera¬ 
ture  functions  representative  of  phenomena  observed  in  actual  shock  tunnel 
operations.  The  mathematical  formulation  used  in  the  analog  studies  was 
adapted  to  the  development  of  an  economical  and  general  digital  data  reduction 
program.  A  description  is  given  of  the  operation  and  design  details  of  a 
compensated  analog  circuit  which  operates  directly  on  the  electrical  output 
of  the  thin-film  gage  to  provide  an  output  signal  proportional  to  the  true 
heating  rate. 


mcvrfm  nos  BUMC— MOT  ram 


iMH-vmiinKr 


TABLE  OF  CONTENTS 

I  INTRODUCTION .  I 

n  SURFACE  THERMOMETRY  AND  DATA  REDUCTION  ....  2 

III  ELECTRICAL  AND  THERMAL  CHARACTERISTICS  OF 

THE  THIN-FILM  RESISTANCE  THERMOMETER .  5 

IV  ANALOG  COMPUTER  STUDIES .  6 

1.  OBJECTIVES .  6 

2.  GENERAL  MATHEMATICAL  FORMULATION .  6 

3.  SPECIFIC  DETAILS  OF  ANALOG  SIMULATION  ....  9 

4.  TEST  PROGRAM  AND  RESULTS . 12 

V  DEVELOPMENT  OF  A  DIGITAL  COMPUTER  PROGRAM  .  .  15 

1.  OBJECTIVES . 15 

2.  MATHEMATICAL  BASIS . 15 

3.  RESULTS . 16 

VI  DEVELOPMENT  OF  A  COMPENSATED  ANALOG  CIRCUIT  .  18 

1.  BACKGROUND  AND  OBJECTIVES . 18 

2.  CORRECTION  FUNCTIONS . 18 

3.  CIRCUIT  DESCRIPTION . 19 

VII  SUMMARY  . 22 

VHI  REFERENCES . 23 

APPENDIX  A  SUMMARY  COMMENTS  ON  ANALOG  COMPUTER  RUNS  33 

APPENDIX  B  DIGITAL  COMPUTER  PROGRAM . 86 

APPENDIX  C  DERIVATION  OF  CORRECTION  FUNCTIONS  FOR  THE 

COMPENSATED  ANALOG  CIRCUIT . 96 

APPENDIX  D  COMPUTER  PROGRAM  FOR  DEFINING  COMPONENT 

VALUES  FOR  THE  ANALOG  NETWORK . 102 


SECTION  I 


INTRODUCTION 


The  thin-film  resistance  thermometer,  comprised  of  a  platinum  alloy 
element  on  a  Pyrex  substrate,  has  been  used  extensively  for  measurements 
of  heat  transfer  in  short-duration  hypersonic  test  facilities  because  of  its 
simplicity  of  principle  and  construction  as  well  as  very  rapid  response.  Since 
the  output  signal  from  the  gage  is  a  measure  of  the  substrate  surface  tempera¬ 
ture  rise,  rather  than  the  incident  heat  transfer  rate,  signal  processing  is 
required  to  obtain  the  desired  parameter.  By  application  of  classical  solid 
conduction  theory,  surface  temperature -time  data  may  be  transformed  into 
the  incident  heat  flux  provided  the  electrical  and  the  thermal  properties  of 
the  gage  are  known. 

For  modest  heating  rates  that  result  in  a  substrate  surface  temperature 
increase  of  the  order  of  100  °F  or  less  during  the  test  period,  the  electrical 
characteristic  of  the  thermometer  may  be  taken  as  a  linear  function  of  tempera 
ture  and  the  substrate  thermal  properties  assumed  invariant.  For  these 
particular  circumstances,  both  analog  and  digital  techniques  for  data  reduction 
have  long  been  available  * »  since  a  closed  form  mathematical  solution  to  the 
problem  exists.  On  the  other  hand,  for  substrate  surface  temperature  excur¬ 
sions  up  to  1000  °F,  the  upper  practical  limit  of  the  thin-film  gage,  the  elec¬ 
trical  characteristic  of  the  resistance  thermometer  is  highly  nonlinear  and 
the  temperature  variation  of  the  substrate  thermal  properties  results  in  a 
nonlinear  form  of  the  heat  diffusion  equation  for  which  no  closed  form  mathe¬ 
matical  solution  exists.  Data  reduction  techniques  in  this  regime  (high  heating 
rates)  are  not  as  well  defined. 


This  report  discusses  studies  directed  toward  the  formulation  of  digital 
and  analog  techniques  for  the  conversion  of  thin-film  resistance  thermometer 
data  to  heat  transfer  rate  with  emphasis  on  high  heating  rates  and  nonconstant 
heat  flux  conditions.  Specifically,  the  program  consisted  of  three  distinct 
efforts:  (1)  analog  computer  solutions  of  a  variety  of  temperature  histories 
which  included  the  study  of  such  factors  as  noise,  tunnel  starting  phenomena, 
and  gage  abrasion;  (2)  development  of  a  practical  digital  computer  program 
for  data  reduction  for  nonspecific  temperature  histories  (and  amplitudes);  and, 
finally,  (3)  formulation  of  a  small,  compact  electronic  analog  circuit  suitable 
for  general  data  conversion  applications. 


The  technology  related  to  the  fabrication  and  application  of  the  thin- 
film  heat  transfer  gage  has  been  presented  in  an  earlier  interim  project 
report-*.  Consequently  only  a  brief  mention  is  made  in  the  next  section  of  the 
underlying  principles  of  the  gage.  These  comments  are  then  followed  by  a 
short  discussion  of  the  electrical  and  thermal  characteristics  of  the  gage 
together  with  numerical  data.  The  three  parts  of  the  program  are  then  pre¬ 
sented  and  the  report  is  concluded  with  a  summary. 


Superscripts  denote  references. 
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SECTION  II 


SURFACE  THERMOMETRY  AND  DATA  REDUCTION 


The  thin-film  heat  transfer  gage  is  a  device  that  provides  a  measure¬ 
ment  of  the  surface  temperature  of  a  substrate  which  is  characterized  as  semi 
infinite,  homogeneous,  and  isotropic.  By  application  of  classical  solid  con¬ 
duction  theory  surface  temperature -time  data  may  be  transformed  into 
the  incident  heat  transfer  rate  when  the  substrate  thermal  properties  are 
known.  To  achieve  a  very  rapid  response  of  the  sensor,  and  to  simplify  the 
data  reduction  procedures,  the  thickness  of  the  thin-film  resistance  thermom¬ 
eter  is  made  sufficiently  small  that  its  presence  on  the  substrate  has  a 
negligible  effect  on  the  substrate  surface  temperature.  The  characteristic 
diffusion  time  for  a  metallic  film  about  four  microinches  thick  is  approxi¬ 
mately  a  microsecond;  hence  for  practical  purposes,  the  film  senses  the 
instantaneous  substrate  surface  temperature. 

An  insulator  material,  in  both  a  thermal  and  electrical  sense,  is  used 
as  the  substrate,  functioning  mechanically  as  a  support  for  the  resistance 
film  and  thermally  as  a  heat  sink.  Sensitivity  in  the  measurement  of  very  low 
heat  transfer  rates  is  enhanced  by  the  use  of  a  substrate  with  poor  thermal 
conductivity.  For  practical  shock  tunnel  applications,  #7740  Pyrex  is  used 
widely  as  a  gage  substrate  so  that  even  in  small  thicknesses  (say  0.04  to  0.  06 
inch)  it  represents  a  semi-infinite  heat  sink  because  of  the  small  depth  of  heat 
penetration. 

The  transient  one -dimensional  heat  conduction  problem  for  a  semi¬ 
infinite  slab  of  material  with  heat  flux  normal  to  the  surface  is  illustrated  by 
the  following  sketch. 


Q  =  heat  transfer  rate 
k  =  thermal  conductivity 
/o  -  density 
c  -  specific  heat 
T  =  temperature 
t  =  time 


For  conditions  of  constant  substrate  properties,  the  thermal  diffusion 
equation  is 

ar(x,t )  _  _a_  <?zT{z,t)  (1 

dt  <oc  gx* 
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with  boundary  conditions 


T(x,o)  =  o  x  >  O 

f(o,t)-  -  r  t  >o 

•  J  -r 


T(x,i  )z  O  ,  X  -*  oo  ,  t  >o 

One  form  of  solution  to  Eq,  (1),  the  derivation  of  which  is  given  in 
detail  in  the  literature  on  thin-film  gages  1,  is  given  below 


i t  (:*)’ i  (**)' 


w.,.  1  r  **  rM-ti  r(\)  ,, 
™*nj  — TTTJV. " 


(2) 


where  T(t)  is  the  substrate  surface  temperature  and  \  is  a  dummy  variable 
of  integration.  This  form  of  the  equation  is  readily  adaptable  for  evaluation 
by  a  digital  computer.  For  the  special  situation  when  the  incident  heat  trans 
fer  rate  is  constant,  Eq,  (2)  reduces  to  the  following  form. 


f  (o,  t)  = 


(TTfic-e*)*  T(t) 
2i  T 


(3) 


The  increase  in  surface  temperature  occurs  in  proportion  to  the  square  root 
of  time  and  the  recorded  gage  output  signals  consequently  are  of  parabolic 
form. 


A  parabolic  surface-temperature  function  has  been  assumed  in  formu¬ 
lating  Eq.  (2).  As  the  temperature  function  deviates  from  this  form,  results 
obtained  from  Eq.  (2)  become  increasingly  less  accurate^.  A  perfectly  general 
expression,  which  does  not  involve  initial  assumptions  about  the  form  of  the 
temperature  function,  is  given  in  Ref.  4. 


f  p  T(tJ-  r(tc.,)  1 

t  -  Trt_/ 

. 

<v 

(4) 


In  deriving  Eq.  (4)  however,  the  temperature  function  is  approximated  by  a 
piecewise  linear  function  so  that  the  accuracy  of  the  result  obtained  is  only 
limited  by  the  degree  to  which  the  actual  function  is  approximated  by  this 
means.  Accuracy  is  improved  therefore  by  choosing  n  sufficiently  large  to 
insure  a  satisfactory  approximation  to  the  true  temperature  function.  Equation 
(4)  is  also  easily  programmed  for  digital  computer  calculations. 

On  the  other  hand,  analog  techniques  have  been  developed  for  data 
reduction  by  taking  advantage  of  the  fact  that  differential  equations  which 
describe  the  one-dimensional  heat  flux  into  a  semi-infinite  solid  are  identical 
with  those  that  apply  to  the  current  flow  in  a  distributed  resistance-capacitance 
line.  A  variety  of  simple,  passive  analog  circuits,  which  operate  in  real 
time,  have  been  developed  for  converting  thin-film  heat  transfer  gage  signals 
directly  into  signals  proportional  to  heat  flux.  One  particular  type,  described 
in  Ref.  2,  has  been  used  extensively  at  CAL.  Except  for  the  need  to  match 
the  frequency  range  of  the  analog  circuit  to  the  Fourier  frequency  components 
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of  the  temperature  function,  there  are  no  assumptions  or  restrictions  concern¬ 
ing  the  nature  of  the  temperature  function. 

For  test  conditions  at  which  the  substrate  thermal  properties  cannot  be 
taken  as  constant,  Eq.  (1)  must  be  recast  in  the  following  form 

-  h  \4(t)£]  (5> 

which  has  no  closed  form  analytical  solution.  Numerical  techniques  may  be 
used  to  evaluate  Eq.  (5)  in  terms  of  the  incident  heat  flux  q(t)  which  produces 
the  measured  surface  temperature  T(t).  The  specific  method  which  was  used 
on  this  project  to  obtain  both  the  analog  computer  and  the  digital  computer 
solutions  to  Eq.  (5)  is  described  in  detail  in  sections  IV  and  V  of  this  report. 
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Hartunian  and  Varwig  have  solved  Eq.  (5)  by  making  use  of  a  trans¬ 
formed  variable,  0  =  /  kdT  and  linearized  forms  for  the  transformed 

Jo 

variable  and  the  thermal  diffusivity.  In  essence,  a  transformed  surface 
temperature  is  obtained  which  when  used  with  the  linear  digital  and  analog 
data  reduction  techniques  described  above  will  yield  the  true  incident  heat 
flux.  A  detailed  account  of  this  procedure,  which  was  used  in  computing 
design  information  for  the  electronic  analog  circuit,  is  included  in  Appendix  C 
of  this  report. 

Finally,  note  should  be  made  of  the  work  of  Schmitz^  in  the  development 
of  a  nonlinear  analog  network  to  directly  compensate  for  the  changes  in  sub¬ 
strate  properties.  The  circuit  proved  to  be  too  complex  for  practical  use  and 
was  further  restricted  by  an  ability  to  compensate  for  changes  in  but  one  of 
the  substrate  parameters  (specific  heat). 

The  purpose  of  the  preceding  discussion  has  been  to  describe,  in  general 
terms,  the  various  methods  used  in  the  computation  of  heat  flux  from  thin- 
film  heat  transfer  gage  measurements  of  surface  temperatures.  It  is  apparent 
that  preceding  this  investigation  relatively  little  emphasis  has  been  placed  on 
solutions  to  the  nonlinear  heat  diffusion  equation  (Eq.  5).  This  consideration 
led  to  the  present  study  of  economical  methods  of  analog  and  digital  data 
reduction  for  the  nonlinear  heat  conduction  situation. 
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SECTION  III 


ELECTRICAL  AND  THERMAL  CHARACTERISTICS  OF  THE 
THIN-FILM  RESISTANCE  THERMOMETER 


Two  distinct  operations  must  be  performed  on  the  output  signal  of  the 
thin-film  resistance  thermometer;  first,  a  correction  must  be  made  for  the 
nonlinear  res istance -temperature  function  of  the  resistance  element  (to  convert 
the  measured  surface  temperature  to  a  true  surface  temperature)  and  second, 
the  variation  of  the  gage  thermal  parameters  with  temperature  must  be  accom¬ 
modated  in  the  data  analysis. 

Measurements  of  the  electrical  and  thermal  characteristics  of  the  thin- 
film  resistance  thermometer  in  the  range  from  room  temperature  of  1000  °F 
have  been  made  at  CAL^.  These  tests  were  performed  on  gages  which  used 
Hanovia  Bright  Platinum  05 -X  solution  for  the  resistance  thermometer  and 
#7740  Pyrex  as  the  substrate. 

The  resistance-temperature  characteristic  of  the  film  is  the  quadratic 
function 


~  -  (T  -TO)  -  A  2“?5  x  f  o'  *  (  t-  70)  ^  -  T  (6) 

m7  0 

AR  =  change  in  gage  resistance 

K70  =  slope  of  the  resistance  function  (dR/dT)  at  70  °F 

T  =  actual  temperature 

Tm  =  measured  gage  temperature  change* 

Table  I  provides  a  quantitative  tabulation  of  the  magnitude  of  the  correction 
required  to  the  temperature  change  as  measured  by  the  gage.  Corrections  are 
of  the  order  of  +15%  at  1000  °F. 


The  temperature  variation  of  the  substrate  thermal  properties  was 
determined  by  an  electrical  pulsing  technique®  that  provided  a  constant  flux 
into  the  resistance  thermometer.  Using  Eq.  (3),  the  value  of  (  p  ck)  as  a 
function  of  temperature  was  obtained.  However,  as  shown  in  Eq.  (5),  the 
parameters  required  are  p  c(T)  and  k(T)  rather  than  (/©ck)  (T).  Handbook 
data  for  p  c(T)  for  #7740  Pyrex  were  used  and  k(T)  was  computed  from  the 
experimental  data**.  Table  II  lists  the  values  for  p  c  and  k  that  were  exclu¬ 
sively  employed  in  all  of  the  studies  included  in  this  report. 


Measured  gage  temperature  change  is  that  obtained  by  assuming  a  linear 
resistance  change  with  temperature  and  using  K-,q  as  the  slope  of  the  function. 

Whereas  published  data  for  p  c(T)  are  in  excellent  agreement,  wide  disagree¬ 
ment  exists  on  measurements  of  k(T)  (see  Ref.  9,  for  example). 
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SECTION  IV 


ANALOG  COMPUTER  STUDIES 


1.  OBJECTIVES 

Analog  computers  are  uniquely  appropriate  to  the  solution  of  partial 
differential  equations  and  offer  advantages  of  economy  and  flexibility  especially 
when  the  effects  of  variations  in  parameters  are  to  be  explored.  Accordingly, 
a  limited  analog  computer  study  was  conducted  at  CAL  several  years  ago  using 
the  computing  equipment  available  at  that  time*.  The  purpose  of  that  study 
was  to  ascertain  the  magnitudes  of  the  corrections  required  to  measure  heat¬ 
ing  rates  as  the  result  of  substrate  property  changes  with  temperature. 
Accuracy  limitations  of  the  analog  equipment  precluded  any  extensive 
investigations. 

With  improved  analog  computer  hardware  now  available  at  CAL,  it  was 
believed  the  project  objectives  aimed  at  developing  performance  criteria  for 
computer  programs  and  analog  network  designs  for  conditions  of  nonconstant 
heating  rates  could  be  furthered  by  solving  the  nonlinear,  one -dimens ional 
heat  diffusion  equation  on  these  facilities.  A  variety  of  temperature  histories 
would  be  used  as  inputs  to  which  would  be  added  such  modifying  factors  as 
tunnel  starting  phenomena,  random  noise,  and  sinusoidal  modulation.  In  all 
cases  the  interest  was  in  the  influence  of  those  disturbances  on  the  heat  flux 
history.  Input  variables  were  limited  to  a  real  life  test  time  of  0.  010  second 
and  a  maximum  temperature  change  of  1000  °F. 

2.  GENERAL  MATHEMATICAL  FORMULATION 


The  nonlinear  diffusion  equation  for  one -dimens ional  heat  conduction  in 
a  homogeneous  and  isotropic  slab  of  material  of  thickness  Ji  is  given  below. 


3T  _  _± _ 

dt  fiC(T) 


(7) 


T  =  temperature 

t  =  time 

x  =  depth  into  slab  (measured  normal  to  surface) 

p  -  density  of  slab 

c  =  specific  heat  of  slab 

k  =  thermal  conductivity  of  slab 

Assumed  boundary  conditions  are  as  follows: 

T  =  T  (x,t) 

T  (0 ,  t)  =  T  (t)  (forcing  function) 

T  (0,0)  =  70°F  (assumed  initial 
T  (i,  t)  =  70°F  temperature) 


The  results  of  this  effort  have  not  been  published. 
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The  last  relation  above  indicates  that  the  slab  is  sufficiently  thick  (semi¬ 
infinite)  that  the  rear  surface  temperature  remains  constant  for  the  duration 
of  the  test  event. 

For  computational  purposes,  Eq.  (7)  is  recast  into  a  different  but 
equivalent  form.  If  we  let  q  represent  the  heating  rate  (in  general  q  =  q  (x,  t)  ), 
then  by  definition 

r  ~  ~  (8) 

$ 

By  substitution,  Eq.  (9)  is  obtained  . 

_  _l  a  /Qv 

c?X  fiC  d% 


This  partial  differential  equation  may  be  approximated  with  a  finite  set 
of  ordinary  differential  equations  developed  by  quantizing  the  space  variable, 
x,  into  n  +  2  values:  xQ,  xj,  ...  xn,  xn  +  p  The  front  surface  of  the  slab 
would  be  denoted  by  x0  and  the  rear  surface  by  xn  +  \  while  xj,  X2,  .  .  .  xn 
would  correspond  to  spatial  locations  within  the  slab  at  which  T  and  q  are  not 
known  (not  specified  by  boundary  conditions). 


To  illustrate  the  concept  used  in  formulating  the  set  of  approximate 
equations,  consider  first  the  case  of  using  a  uniform  spacing  in  the  slab  depth 
x.  Let  x  take  the  values  o,  h,  2h,  3h,  ...  (n  +  l)h  with  n  odd  and,  as  a  matter 
of  choice,  let  the  value  of  q  be  calculated  at  the  odd  multiples  of  h  and  T  values 
be  calculated  at  even  values  of  h.  The  figure  below  illustrates  this  concept. 


X 


Tq  and  Tn  +  ^  are  specified  boundary  conditions 
Thus  the  equation 


may  be  approximated  as 


7z  -  To 
2  h 


where  =  k  (T^)  . 


This  substitution  has  made  use  of  the  linear  form  of  Eq.  (7). 

k^  must  be  evaluated  in  terms  of  Tq  and  T2  since  T  values  are  not  defined 
for  an  odd  subscript.  In  actual  practice  k.  was  taken  as  the  value  of  k  at 

Jo  . 
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The  companion  equation  for  T 


[^t\  =  ['k  W] 


may  be  approximated  analogously  as 


IZk  = .  JL  /JiJl.) 

dt  (  */i  / 


where  =  T^  (t)  =  T  (x2»t) 

'oc2  =  'oc(t»*=*2  =  <T2> 

By  this  procedure,  the  following  set  of  equations  is  generated. 


in  = 


+  1  ~Tti-  / 


,  *  =  3 ,  5  •  •  . 


(10a) 


$  »  +  /  ~  9  n  -  f 


,  ^  *2,^/ 6 


(10b) 


Thus  with  T0  and  Tn  +  j  known  and  T2,  T4  .  .  .  Tn  _  j  the  unknown  dependent 
variables  (with  the  "q's"  algebraically  related),  Eq.  (10a)  and  (10b)  represent 

a  set  of  -7>--  —  simultaneous  differential  equations.  These  equations  are 

quadratically  correct  in  the  sense  that  if  (for  a  fixed  t)  T  (x,  t)  varies  at  most 
quadratically  in  x,  Eq.  (10b)  is  exact.  Equation  (10a)  and  (10b)  may  be  regarded 
as  the  results  obtained  when  correction  terms  proportional  to  h^,  h  ,  etc.  are 
dropped  (the  validity  of  this  procedure  is  dependent  on  h  being  made  sufficiently 
small). 

The  result  desired  from  the  solution  of  Eq.  (10a)  and  (10b)  is  the  surface 
heat  transfer  rate,  qQ  =  q  (0,t)  whereas  the  direct  results  obtained  provide 
qj,  q3»  etc.  at  interior  planes.  If  h  is  small,  the  assumption  may  be  made 
that  qj  =  q  .  An  appreciable  increase  in  accuracy  is  easily  achieved  however 
by  appropriate  extrapolation  to  the  surface  heating  rate. 

An  analysis  of  the  errors  can  be  made  by  assuming  qj  =  q0.  Taking 
the  linear  diffusion  equation,  the  following  well-known  relation  is  obtained  as 
a  solution. 


T(d,i) 


Kt  r 
rt  l- 


i~Ki 


etfc 


*=■] 

let  J 


2//FF 


thermal  diffusivity 


k//oc 
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Differentiating  Eq.  (11)  yields  the  following  relation  for  the  temperature 
gradient  and  hence  heating  rate^. 


j. 

*  dx 


*<*•*)• 


(12) 


At  a  fixed  t,  Eq.  (12)  is  nearly  lineal  in  x  for  small  values  of  x. 


=  fa  ['  ~ 


X 

VTTt 


(13) 


Thus  a  linear  extrapolation  formula  is  reasonable  for  the  case  of  the  linear 
diffusion  equation  and  by  inference  for  the  nonlinear  case  as  well  since  the 
thermal  substrate  properties  are  smoothly  changing  functions  of  T.  A  much 
better  estimate  than  q^  by  itself  is  obtained  by  use  of  Eq.  (14). 


a  s*lLJs 

To  2 


(14) 


While  Eq.  (10a)  and  (10b)  are  conceptually  simple,  the  use  of  uniform 
spacing  in  x  over  the  entire  thickness  of  the  slab  is  not  practical.  To  achieve 
accuracy,  the  spacing  h  would  have  to  be  very  small  so  that  the  temperature 
variation  over  h  is  also  small.  Since  most  of  the  temperature  variation  occurs 
very  near  the  front  surface  (in  materials  of  relatively  low  thermal  diffusivity 
like  Pyrex),  a  value  of  h  of  the  order  of  5  microinches  would  be  implied. 

Thus,  for  slabs  fulfilling  semi- infinite  depth  requirements  in  shock  tunnel 
applications,  several  thousand  differential  equations  would  be  required.  Of 
this  large  total,  relatively  few  would  have  physical  significance. 


Two  equally  logical  modifications  are  available  to  avoid  the  problems 
associated  with  a  uniform  spacing  in  x.  One  approach  is  to  use  nonuniform 
x  spacing  (requiring  replacement  of  Eq.  (10a)  and  (10b)  with  more  complicated 
expressions)  while  the  other  involves  a  change  of  variable,  e.  g.  ,  x  =  7  z^  or 
x  =  7  z^,  etc.  with  uniform  spacing  in  the  new  variable.  The  latter  approach 
was  considered  most  expedient  and  therefore  adopted. 


3. 


SPECIFIC  DETAILS  OF  ANALOG  SIMULATION 


The  earlier  (unreported)  analog  computer  studies  had  made  use  of  the 
transformation  x  =  y  z^  (  y  =  constant)  which  was  chosen  somewhat  arbitrarily. 
A  review  was  made  accordingly  at  this  time  to  determine  the  appropriateness 
of  this  transformation. 


As  stated  previously,  near  the  front  surface  of  the  slab  (which  is  the 
region  of  most  significance),  q  is  nearly  linear  in  x  so  that  it  would  be  cubic 
in  z  (q  aQ  +  a^x  ~  aQ  +  a^z^).  The  simpler  forms  of  difference  equations, 
analogous  to  Eq.  (10a),  however  are  not  accurate  for  cubic  terms  but  only 
accurate  for  terms  through  quadratic.  For  this  reason  the  cubic  transformation 
is  not  mathematically  sound  and,  in  effect,  over-emphasizes  the  distortion  in 
the  space  coordinate.  This  effect  has  also  more  than  mathematical  implica¬ 
tions.  Rather  severe  scale  factor  and  accuracy  problems  are  encountered 
in  analog  computer  implementation.  In  effect,  very  small  voltages  are  formed 
by  the  difference  in  two  large  (nearly  equal)  voltages.  These  low  signal  levels 
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also  result  in  unfavorable  signal-to-noise  ratios.  Consequently,  a  trans¬ 
formation  nearly  linear  for  small  x,  and  at  most  quadratic,  is  desirable. 

The  x  =  ?■  z6  (  y  =  2,5  x  10”  )  transformation  was  chosen  for  analog 
computation  and  very  satisfactory  locations  of  internal  planes  were  thereby 
achieved.  High  resolution  near  the  surface  was  realized  simultaneously  with 
adequate  coverage  of  the  deeper  internal  regions.  Furthermore,  the  trans¬ 
formation  satisfies  the  "quadratic"  requirement  of  simple  difference  approxi¬ 
mations  to  derivatives.  Locations  and  designations  of  the  planes  defined  by 
this  transformation  are  tabulated  below. 

Temperature  and  Heating  Rate  Plane 
Locations  in  Actual  and  Transformed  Spatial  Coordinates 


microfeet  microfeet 


Plane 

T  =  temperature 
4  =  heating  rate 


2.  5 

1 

.  0 
«1 

10.  0 

2 

T 

A  ? 

22.  5 

3 

9* 

L 

40.  0 

4 

J 

T4 

62.  5 

5 

*5 

90.  0 

6 

T6 

122.  5 

7 

*7 

160.0 

8 

T8 

202.  5 

9 

9q 

250.  0 

10 

7 

T10 

360.0 

12 

*12 

490.  0 

14 

T14 

902.  5 

19 

*19 

1440.  0 

24 

T24 

Regardless  of  the  change  in  variable,  it  is  recognized  that  extrapolation 
to  heat  transfer  rate  at  the  surface,  q0,  must  be  based  on  linearity  in  x,  not 
linearity  in  the  artificial  variable  z.  Thus 

x  =  2.  5  x  10  *  z^  feet  =  2,5  z^  microfeet 
so  that  computer  implementation  produced  voltages  corresponding  to 
=  q  (evaluated  at  z  =  1)  =  q  (evaluated  at  x  =  2.  5) 


q^  s  q  (evaluated  at  z  =  3)  =  q  (evaluated  at  x  =  22.  5). 
Thus  to  obtain  qQ  (x  =  0,  z  =  0),  the  correct  relationship  is 


Note  that  the  value  of  qj  is  evaluated  at  a  point  physically  9  times  closer  to 
the  surface  than  is  q^,  hence  the  9:1  relative  weighting. 
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A  listing  of  the  set  of  linear  difference  equations  and  the  ordinary  differential 
equations  actually  used  is  given  below.  A  tabulation  of  the  values  used  for 
k(T)  and  ^oc(T)  is  shown  in  Table  II. 


Difference  Equations  for  q 


.  * 


i<  = 

4s  * 
U  s 
h  = 


= 


*/  = 


'/« lift 

■  4s  ^ 

.  *1 

TZ-To 

zrh 

2 

-  *3 

7V-T2 

2 

-  *5 

r6  -  T4 

2S'?S 

2 

-*7 

T*-T<. 

2??7 

2 

-*9 

T/d  * 

2 

-*I2 

Tl+'Tlo 

f 

-■*/9 

T2A  " 

*»7/? 

to 

A(tz)+A(t0)  ' 

o  j  ** 

/ 


10  x  10 

_ / 

30  x  10 
I 


-<• 


*,  (TZ  -  To ) 

H  *3  (T+  ~  V 
h  <J,  -  T.) 


70  x  10  6  7  y 

/ 


<70  x  io 
I 


24-Ox  10 
I 


izr  *?  <Tio  -  V 
TZ,  ft/2  Tto) 


150  x  10 

ACT+)  +A(TZ) 


-(,  A?  ^T2f  T/V-) 


* 
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See  Eq.  (10a). 
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Differential  Equations  for  T 
Tq:  the  input  function 
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i 
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dr 
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3 
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/ 
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dT 

C!H- 

/O 
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T 2^  =  70  by  hypothesis 

c  =  fic  _3 

t  =  OCX  ;  0C  =  10  (change  from  real  time  to  computer  time) 

4.  TEST  PROGRAM  AND  RESULTS 

A  total  of  46  test  runs  was  made  during  the  program.  In  general,  the 
input  temperature  rise**  was  permitted  to  reach  a  maximum  value  of  1000 °F 
in  10  msec  to  explore  the  full  effects  of  the  variation  in  substrate  thermal 


See  Eq.  (10b). 

It  was  not  considered  important  to  the  results  of  this  study  that  the  correc¬ 
tion  of  nonlinearity  of  the  resistance  characteristic  of  the  thermometer  be 
included  in  the  program.  Accordingly,  the  input  temperature  to  the  analog 
computer  was  taken  as  the  actual  (or  true)  surface  temperature. 
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properties.  Except  in  special  situations  where  the  ease  of  discerning  the 
effects  of  changes  in  the  input  function  would  have  been  compromised,  each 
set  of  test  conditions  was  used  to  generate  two  heat  transfer  rate  histories: 
one  with  the  substrate  properties  constant  and  the  other  with  the  properties 
variable  with  temperature. 

Data  were  recorded  on  a  single  channel  x-y  plotter  and  comprised  the 
temperature  input  signal  and  the  corresponding  heating  rate.  The  resultant 
plots  are  shown  labeled  as  Run  #1  through  Run  #46  and  are  included  as  part 
of  Appendix  A.  Each  record  includes  a  brief  tabulation  of  the  test  conditions 
imposed.  To  an  extent,  these  records  are  self-explanatory;  however,  pertin- 
rrf  details,  comments,  and  observations  applicable  to  each  record  are  included 
Ani.i  ndix  A. 

These  records  quantitatively  demonstrate  the  effects  of  such  factors  as 
variable  substrate  properties,  electronic  filtering,  noise  (random  and  sinu¬ 
soidal),  simulated  tunnel  starting  phenomena,  and  gage  abrasion  on  the  derived 
heat  transfer  rate.  Basically  the  fundamental  input  has  been  a  parabolic 
function  with  lesser  emphasis  given  to  ramp  and  step-type  inputs  as  well. 

For  large  incremental  temperature  changes,  neglect  of  the  temperature 
dependent  thermal  properties  of  the  substrate  (#7740  Pyrex)  results  in  heat 
transfer  rates  that  are  too  low  by  as  much  as  15%.  Tunnel  starting  phenomena 
that  are  of  an  impulsive  nature  (very  short  in  time  extent)  cause  only  a  tran¬ 
sient  disturbance  in  the  heating  rate  and  the  disturbance  dies  out  very  rapidly. 
Other  phenomena,  as  exemplified  by  gage  abrasion,  that  result  in  step-type 
inputs  (that  is  of  long  time  duration)*  also  result  in  a  large  transient  distur¬ 
bance  in  heating  rate.  This  type  disturbance  however  decays  as  t “1/2  and 
hence  has  a  long  persistence  and  results  in  erroneous  data  if  the  q  attributable 
to  this  type  input  is  indeed  spurious  and  not  desired  in  the  final  result.  Elec¬ 
tronic  filtering  of  the  temperature  signal,  often  employed  in  practice,  is  very 
effective  in  reducing  excessive  amplitude  excursions  in  q  resulting  from  noise 
and  step  disturbances  in  the  temperature  signal.  Filtering  however  is  not  at 
all  effective  in  suppressing  the  long  term  (t"^'^)  contributions  of  step  inputs  to 
the  net  q.  The  relative  effect  of  step  inputs  on  q  depends  on  the  relative 
amplitudes  of  the  fundamental  signal  (say  a  parabola)  and  the  spurious  (assumed) 
step.  For  approximately  equal  amplitudes,  the  q  component  corresponding  to 
the  step  signal  maybe  as  much  as  a  factor  of  2  larger  than  the  component  of  q 
corresponding  to  the  parabola  (see  Run  #34). 

The  results  of  these  analog  studies  highlight  the  dichotomy  of  the  problem 
of  heat  transfer  data  reduction.  First  is  the  necessity  of  correctly  trans¬ 
lating  the  temperature  data  into  heat  flux  data  and  second  is  the  extraction  of 
the  desired  component  of  the  resultant  heat  flux  which  may  be  a  composite  of 
two  or  more  components. 


Despite  the  fact  that  the  impact  of  a  particle  on  a  gage  is  a  transient  of  very 
brief  duration  in  a  physical  sense,  the  electrical  effect  is  of  long  duration 
since  a  step-type  increase  in  gage  signal  occurs  and  persists  throughout  the 
test  event. 
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In  the  situation  where  gage  abrasion  exists,  it  is  obvious  that  the  gage 
output  signal  does  not  represent  the  surface  temperature  history  since  the 
signal,  in  part,  represents  resistance  changes  caused  by  mechanical  effects 
rather  than  thermal  effects.  Thus  erroneous  heat  flux  data  may  be  obtained 
despite  the  fact  that  the  data  reduction  procedures  have  been  properly  executed. 
As  illustrated  by  the  data  of  Appendix  A,  the  t"1'^  decay  characteristic  of  the 
heat  flux  corresponding  to  a  step  temperature  change  ensures  the  fact  that  this 
spurious  signal  will  persist  throughout  the  length  of  the  typical  shock  tunnel 
test.  Where  serious  abrasion  exists,  the  only  recourse  is  to  record  the  gage 
output  signal  and  to  reconstruct  the  actual  temperature  signal  by  subtracting 
from  the  gage  signal  the  components  attributable  to  abrasion  (these  instances 
are  usually  manifested  as  sharp,  vertical  discontinuities  in  the  signal).  The 
reconstructed  temperature  signal  is  then  used  in  computing  the  heat  transfer 
rate. 


Whereas  the  effects  of  abrasion  are  well  understood,  the  extent  to  which 
heat  flux  data,  as  computed  for  the  steady  flow  portion  of  the  test,  reflect 
effects  possibly  associated  with  tunnel  starting  phenomena  is  not  well  known. 

If  these  effects  are  of  an  impulsive  character,  their  disturbance  of  the  heat 
flux  signal  will  be  a  short-lived  transient  of  small  concern.  On  the  other  hand, 
should  the  flow  be  such  to  create  a  sharp  increase  in  gage  temperature  which 
then  is  subsequently  further  increased  relative  to  that  level  in  response  to  the 
subsequent  steady  flow  heat  flux,  then  heat  flux  record  interpretation  would  be 
a  problem  (see  Run  #15,  Appendix  A).  It  would  be  a  problem  in  the  sense  that 
although  the  computed  heat  flux  would  in  fact  be  that  actually  experienced  by 
the  gage  (or  model),  only  that  part  of  the  heat  flux  (and  surface  temperature 
rise)  resulting  from  the  steady  flow  portion  of  the  test  is  really  being  sought. 

The  seriousness  of  the  problems  cited  above  is  directly  related  to  the 
relative  magnitudes  of  the  spurious  effects  compared  with  the  steady  state 
phenomena  being  measured.  There  is  a  body  of  experimental  data  which 
suggests  that  for  some  types  of  test  conditions,  measured  heat  flux  data  may 
be  drastically  affected  by  some  unexplained  circumstance.  In  such  situations 
heating  rates  are  much  higher  than  expected  and  have  a  characteristic  time 
function  that  suggests  that  tunnel  starting  conditions  are  involved.  Additional 
studies  are  needed  to  resolve  these  particular  problems. 
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SECTION  V 


DEVELOPMENT  OF  A  DIGITAL  COMPUTER  PROGRAM 
1.  OBJECTIVES 

The  object  of  this  portion  of  the  project  was  to  develop  an  economically 
practical  digital  computer  program  for  solving  the  nonlinear  heat  conduction 
problem.  Specifically,  the  program  was  to  be  designed  for  the  routine  con¬ 
version  of  the  data  measured  with  thin-film  resistance  thermometers  into 
heating  rates.  A  perfectly  flexible  program  was  desired  that  would  combine 
a  short  running  time,  freedom  from  any  restrictions  on  the  functional  form  of 
the  input  data  (gage  surface  temperature  histories),  and  a  high  accuracy. 

Output  data  from  this  program  were  also  to  serve  as  the  standard  for  evaluation 
of  the  operational  performance  of  the  compensated  analog  network  (see  Section 
VI). 


A  further  requirement  of  the  project  was  the  placing  of  the  digital  data 
reduction  program  on  an  operational  basis  on  the  IBM  7094  system  at  WPAFB. 
The  CAL  Computing  Center  currently  employs  an  IBM  System  360/65  Model 
H  and  the  program  was  initially  written  for,  and  evaluated  on,  this  particular 
system. 

2.  MATHEMATICAL  BASIS 

Based  on  the  successful  experiences  with  the  methods  employed  in  solving 
the  nonlinear  diffusion  equation  on  the  analog  computer,  it  was  decided  to  adapt 
the  identical  methematical  approach  to  the  development  of  a  digital  program 
for  data  reduction.  Advantage  was  taken  of  certain  unique  capabilities  of  a 
digital  computer  to  improve  the  accuracy  of  the  computed  data.  Because  of 
the  fact  that  an  appreciable  amount  of  the  mathematics  parallels  that  of  Section 
IV-2  and  since  the  specific  details  are  included  in  Appendix  B,  this  discussion 
will  be  very  general  and  will  emphasize  the  unique  aspects  of  the  digital  formu¬ 
lation  as  contrasted  to  the  analog  formulation. 

As  in  the  analog  computer  application,  the  substrate  was  divided  into  a 
series  of  slabs  by  planes  oriented  parallel  to  the  surface  plane.  If  we  desig¬ 
nate  the  surface  as  the  zero  plane,  then  again  the  even-numbered  planes 
(including  zero)  represent  the  temperature  planes  while  the  intermediate  (odd- 
numbered  planes)  represent  the  heat  flow  planes.  A  set  of  equations  may  then 
be  written  for  each  of  these  planes;  heat  flow  being  expressed  in  terms  of  a 
first  order  difference  formula  in  temperature  (Eq.  10a)  and  the  time  derivative 
of  temperature  as  a  first  order  difference  formula  in  heat  flow  (Eq.  10b).  A 
simple  Runge-Kutta  integration  scheme  was  adopted  to  solve  this  set  of 
equations  since  in  the  present  circumstance  the  time  variable  must  be  quantized 
in  addition  to  the  spatial  coordinate  x  (depth  into  the  substrate). 

Whereas  the  required  nonlinear  spacing  of  the  planes  followed  a  quadratic 
function  (x  =  az^)  in  the  analog  studies,  considerably  more  freedom  in  select¬ 
ing  the  transformation  function  was  new  possible.  Accordingly  a  polynomial 
transformation  (x  =  a^z  -f  a^z^  +  a^z^)  was  employed  together  with  a  larger 
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total  number  of  planes  (twenty  temperature  planes  compared  with  eight  in  the 
analog  program).  The  accuracy  of  the  program  improves  with  an  increasing 
number  of  planes.  Attention  is  called  to  the  selection  of  the  exact  form  of  the 
transformation  which  is  purely  arbitrary.  Despite  the  arbitrariness  entering 
this  selection,  the  essential  facts  are  that  the  desired  linearity  of  spacing 
near  the  surface  (see  Section  IV-3)  is  achieved  and  yet  the  entire  significant 
depth  of  the  slab  is  covered  with  a  reasonably  small  number  of  planes. 

A  change  of  time  variable  was  also  introduced  to  avoid  a  singularity 
near  the  substrate  surface.  This  singularity  is  associated  with  the  necessity 
to  evaluate  the  time  derivative  of  temperature  (Eq.  10b)  which  goes  to  infinity 
at  time  equal  zero  for  a  parabolic  temperature  history.  A  transformation  of 
the  form  t  =  ?  (1  -  e.~az  )  was  made  so  that  t  a.  fit*  for  small  values  of 
fiX  and  t  ^  T  for  large  /3  Z 

The  digital  program  also  includes  a  routine  for  correcting  the  gage 
measured  temperature  to  actual  surface  temperature  (Eq.  6). 

A  more  detailed  description  of  the  digital  program  formulation  is  given 
in  Appendix  B  together  with  a  listing  of  the  program  in  Fortran  IV  language. 

3.  RESULTS 

Proper  operation  of  the  digital  computer  program  was  ascertained  by 
the  use  of  input  data  functions  whose  solutions  could  be  hand  calculated  with 
ease  and  accuracy.  Thus  surface  temperature  histories  of  parabolic  form 
were  employed  since  these  correspond  to  a  constant  level  of  incident  heat 
flux  if  the  gage  properties  are  held  constant.  Adequacy  of  the  planar  spacing 
and  the  stability  of  the  computer  program  is  easily  verified  by  this  process. 

In  computation  techniques  of  this  type  where  time  and  spatial  parameters  are 
quantized,  instability  may  result  when  the  relative  sizes  of  the  increments 
selected  do  not  satisfy  stability  criteria. 

When  program  operation  had  been  properly  verified,  a  series  of  para¬ 
bolic  temperature  functions,  of  various  amplitudes,  was  used  as  a  trial 
problem.  Solutions  were  obtained  both  for  the  case  of  constant  gage  properties 
and  for  the  case  of  temperature  dependent  gage  properties.  Some  of  these 
results  are  summarized  in  Table  III  which  tabulates  the  surface  heat  flux  data 
for  the  two  cases  as  a  function  of  time.  As  in  the  analog  studies,  it  is  obvious 
that  most  of  the  effects  of  changing  properties  are  experienced  in  the  first 
few  milliseconds  since,  for  a  parabolic  function,  the  most  rapid  change  in 
temperature  also  occurs  in  that  time  interval.  That  the  corrections  may  be 
of  significant  proportions  is  illustrated  below  (data  taken  from  Table  III). 

Uncorrected  Corrected 
Time  Heat  Flux  Heat  Flux 

sec.  Btu/ft^sec  Btu/ft^sec  %  Correction 


.  005 

200.  0 

222. 07 

11.0 

.  010 

200.  0 

231. 08 

15.  5 

.  005 

400.  0 

484.  57 

21.  1 

.010 

400.  0 

512.  74 

28.  2 

.  005 

600.  0 

777.81 

29.  6 

.010 

600.  0 

847. 79 

41.  3 
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The  data  reduction  program,  as  presently  constituted  (see  listing  in 
Appendix  B),  is  designed  to  provide  a  high  accuracy  of  computation  by  virtue 
of  the  large  number  of  temperature  and  heat  flow  planes  that  are  employed. 

As  demonstrated  by  the  analog  computer  studies,  wherein  less  than  one  half 
of  this  number  of  planes  were  successfully  employed,  a  smaller  number  of 
planes  could  be  used  if  lesser  accuracy  were  desired  in  exchange  for  a  shorter 
computing  time. 

The  computer  program  also  has  application  to  the  general  study  of  non¬ 
linear  heat  conduction  in  materials  since  it  evaluates  the  temperature  history 
and  the  heat  flux  history  at  many  planes  within  the  body  of  the  material.  In 
its  present  form  however,  the  only  data  printed  out  are:  surface  temperature 
in  °F  (the  input)  and  surface  heat  flux  in  Btu/ft^sec  (the  output)  for  a  period 
of  10  msec  at  0.  05  msec  intervals. 

The  computer  program,  in  card  deck  form,  has  been  delivered  to 
AFAPL  and  some  identical  preliminary  trial  problems  have  been  run  on  both 
the  IBM  7094  system  at  WPAFB  and  the  CAL  IBM  360  Model  65  system  at 
CAL.  A  very  small  discrepancy  in  the  first  and  second  decimal  positions  of 
the  computed  heating  rate  data  has  been  noted  between  the  two  computers. 

At  the  writing  of  this  report  the  basis  for  this  discrepancy  has  not  been 
positively  identified  although  it  is  thought  to  be  associated  with  lesser  precision 
of  the  CAL  computer  which  uses  fewer  bits  per  word. 
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SECTION  VI 

DEVELOPMENT  OF  A  COMPENSATED  ANALOG  CIRCUIT 


1.  BACKGROUND  AND  OBJECTIVES 

Desirable  in  most  experimental  endeavors,  but  especially  important  and 
critical  in  hypersonic  shock  tunnel  testing,  is  immediate  access  to  the  results 
of  measurements  of  the  test  parameter.  Since  the  output  signal  of  the  thin- 
film  heat  transfer  gage  is  a  surface  temperature,  the  common  practice  in 
earlier  years  involved  manual  tabulation  of  temperature-time  data  from  an 
oscilloscope  photograph  of  the  gage  signal.  These  data  were  subsequently 
transferred  to  punched  cards  and  processed  on  a  digital  computer  to  obtain 
the  heat  transfer  rate.  This  procedure  was  very  time  consuming  and  uneconom¬ 
ical.  Despite  the  current  availability  of  instrumentation  capable  of  directly 
converting  the  gage  output  signal  into  a  computer  compatible  format,  the  time 
interval  between  test  and  the  availability  of  the  test  results  is  still  often 
unsatisfactory. 

Because  they  provide  a  direct,  on-the-spot  conversion  of  the  heat  transfer 
gage  output  signal  into  a  signal  proportional  to  heat  flux,  analog  networks  have 
found  a  wide  acceptance  in  shock  tunnel  testing.  The  principal  drawback  of  the 
simple,  passive  analog  network  is  that  it  neither  can  compensate  for  the  non¬ 
linearity  of  the  gage  resistance-temperature  function  nor  the  effects  of  changing 
substrate  thermal  properties.  Despite  development  of  several  compensated 
electronic  analog  circuits”*  little  practical  use  appears  to  have  resulted. 

Practice  at  CAL  has  been  to  employ  the  simple  linear  networks  and  (where 
appropriate)  to  manually  apply  corrections  to  the  measured  heating  rate.  The 
applicable  corrections,  derived  from  digital  computer  solutions  for  parabolic 
type  temperature  histories,  are  applied  as  functions  of  time  and  measured 
heat  flux. 

In  these  less  than  satisfactory  circumstances,  a  need  would  seem  to 
exist  for  an  economical,  reliable  and  compact  compensated  analog  network 
that  could  be  used  in  each  thin-film  thermometer  channel  to  provide  a  direct 
output  of  the  true  heat  transfer  rate.  Development  of  a  circuit  design  aimed 
at  achieving  these  objectives  was  therefore  undertaken. 

2.  CORRECTION  FUNCTIONS 

As  described  in  Section  III,  two  distinct  and  separate  operations  are 
needed  on  the  output  signal  from  a  thin-film  heat  transfer  gage  if  the  true 
heat  flux  is  to  be  computed.  These  operations  include  the  corrections  for  the 
nonlinear  resistance-temperature  function  of  the  resistance  element  as  well 
as  the  changes  in  substrate  properties. 

The  functional  form  of  the  resistance  nonlinearity  is  given  in  Eq.  (6) 
and,  in  tabular  form,  in  Table  I.  A  mathematical  relationship,  based  on  gage 
output  voltage  and  appropriate  to  circuit  design  considerations,  is  derived  in 
Appendix  C  and  given  below. 

eM  -  0.320  x  10*  -  aJ  OJ+9*  /O9-  o.  077 /  X  IQ*  toy  (15) 
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=  temperature  rise,  °F 

-  gage  voltage  rise 

Uq  -  initial  gage  voltage 

To  derive  a  compensation  function  for  the  changes  in  substrate  properties 
with  temperature,  the  mathematical  techniques  of  Hartunian  and  Varwig^  were 
employed.  A  solution  to  the  nonlinear  diffusion  equation  is  obtained  as  a  trans¬ 
formed  temperature  rise  expressed  in  terms  of  the  true  surface  temperature 
rise.  The  transformed  temperature  represents  the  temperature  that  would 
have  been  experienced  by  the  substrate  had  its  properties  remained  constant. 
Thus  the  availability  of  the  transformed  temperature  permits  the  use  of  the 
simple,  stable,  reliable  passive  analog  network  to  obtain  the  true  heat  transfer 
rate. 


The  expression  for  the  transformed  temperature,  0^,  ,  based  on  the 
thermal  property  data  of  Table  II  is  derived  in  Appendix  C. 

0^  =  0^1  1  +  .0001882  (5/9)  ]  (16) 

Equations  (15)  and  (16)  specify  the  operations  required  on  the  gage  output 
signal  and  form  the  basis  for  the  design  of  the  function  generator  of  the  elec¬ 
tronic  circuit. 

3.  CIRCUIT  DESCRIPTION 

The  nonlinear  operations  (Eqs.  15  and  16)  used  to  convert  the  gage  out¬ 
put  signal  to  an  analog  of  the  transformed  temperature  may  be  combined  into 
a  single  function  and  hence  accommodated  by  a  single  circuit.  A  diode-type 
function  generator  was  selected  for  this  purpose  and  is  shown  as  network  No.  1 
in  Fig.  1  which  is  a  schematic  of  the  entire  circuit.  Figure  2  is  a  simplified 
schematic  diagram  and  will  be  useful  in  graphically  illustrating  the  following 
verbal  description  of  the  circuit  operation. 

In  this  schematic  (Fig.  2)  it  is  presumed  that  the  gage  signal  rises  from 
a  zero  level  and  that  diodes  are  ideal  conductors  in  the  forward  direction  and 
perfect  open  circuits  in  the  reverse  direction.  Voltage  limits  are  also  ignored. 
A  30-milliampere  current  source  produces  a  potential  difference  of  3  volts 
across  a  100-ohm  gage  at  room  temperature.  With  heat  flux  applied  to  the 
gage,  its  resistance  increases  and  the  voltage  at  the  input  of  amplifier  I  sweeps 
positively  reaching  approximately  +3  volts  at  a  gage  temperature  of  1000  °F 
above  room  temperature.  The  amplifier  output,  five  times  as  large  as  the 
input,  sweeps  negatively  by  15  volts  and  is  applied  to  amplifier  II  by  Rj.  As 
this  voltage  becomes  increasingly  negative,  the  diodes  that  are  attached  to  the 
voltage  divider  string  at  the  output  of  amplifier  I  begin  to  conduct  at  intervals 
of  one  volt  and  shunt  R^  with  the  resistors  in  series  with  their  anodes.  Current 
flow  through  these  resistances  into  the  summing  junction  (J2)  of  the  second 
amplifier  increases  in  a  nonlinear  fashion  with  the  output  voltage  of  amplifier  I. 
The  output  of  amplifier  II,  proportional  to  this  current,  will  vary  nonlinearly 
as  shown  in  the  inset  graph  in  Fig.  2. 

Series  resistances  in  the  diode  network  are  selected  to  generate  the 
composite  compensation  function  obtained  by  eliminating  between  Eqs.  (15) 

and  (16)  and  expressing  <f>  in  terms  of  U„  !u„  .  Whereas  it  is  possible  to 

0  00 
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calculate  exactly  the  various  resistances  needed  for  this  purpose,  minor 
variations  in  diode  properties  require  corrections  of  a  few  percent.  There¬ 
fore,  it  is  realistic  to  calculate  only  the  approximate  values  of  the  resistances 
and  adjust  each  resistor,  one  at  a  time,  to  fit  the  desired  function  precisely. 
Adjustment  is  simplified  when  the  diodes  come  into  conduction  at  intervals  of 
more  than  one -half  volt;  however,  the  accuracy  of  curve  fitting  is  poorer  the 
fewer  the  diodes.  In  the  present  application,  a  network  of  six  diodes  at  one- 
half  volt  intervals  proved  less  satisfactory  than  one  with  fifteen  diodes  oper¬ 
ating  at  one  volt  intervals.  This  latter  network,  schematically  shown  in  Fig.  3, 
required  a  larger  voltage  swing  from  the  amplifier. 

By  actual  trial,  only  thirteen  biased  diodes  were  required  to  generate 
the  required  compensation  function  since  two  of  the  diode  series  resistances 
were  determined  to  assume  infinite  values  (i.  e.  ,  the  respective  diodes  need 
to  be  open  circuit).  Bias  voltage  applied  to  the  entire  function  generator  is 
adjustable  by  means  of  the  trimming  potentiometer  shown  in  Fig.  3. 

The  diode  function  generator  provides  the  proper  compensation  to  gage 
measured  temperature  changes,  Tj^*,  up  to  a  maximum  of  1000  °F  which 
temperature  rise  is  taken  to  correspond  to  a  3-volt  signal  output  from  the  heat 
transfer  gage.  As  a  result  of  the  5X  gain  of  amplifier  I,  the  function  generator 
experiences  a  15-volt  signal  in  response  to  a  3-volt  signal  from  the  gage.  In 
the  manner  described  earlier,  the  function  generator  augments  the  applied 
signal  in  a  nonlinear  fashion  with  increasing  signal  amplitude.  The  input- 
output  relationship  is  qualitatively  shown  in  the  inset  graph  on  Fig.  2. 

Calibration  of  the  function  generator  may  be  accomplished  statically  by 
simulating  the  gage  output  signal  )  in  discrete  steps  in  the  range  from  0 
to  3  volts  and  measuring  the  output  signal  ^  .  The  required  functional 
relationship  between  and  is  given  by  Eqs.  (15)  and  ( 1 6) .  A  linear 
time  dependent  signal  (a  ramp)  may  be  used  if  a  dynamic  check  is  preferred. 

In  either  case,  the  Mt  -  function  will  have  the  form  shown  in  the  inset 
graph  of  Fig.  2. 

The  output  voltage  range  of  the  amplifier  should  be  as  large  as  practically 
possible  to  improve  the  accuracy  of  curve  fitting  and  to  improve  the  signal-to- 
noise  ratio.  Analog  computer  amplifiers  normally  possess  a  nearly  sym¬ 
metrical  range  of  output  about  zero  volts.  It  is  therefore  advisable  to  bias  the 
amplifiers  so  that  the  full  voltage  range  can  be  used  with  an  input  signal  that 
is  always  of  one  sign.  A  number  of  resistors  shown  in  Fig.  1,  but  not  in 
Fig.  2,  serve  this  purpose. 

Two  diodes  at  the  output  of  amplifier  I  introduce  a  distortion  into  the 
voltage  applied  to  the  diode  network  which  essentially  cancels  the  nonideal 
logarithmic  characteristic  of  the  network  diodes.  A  direct  feed  from  the  end 
of  the  feedback  resistor  of  amplifier  I  to  the  input  of  amplifier  II  assures  that 
this  distortion  is  not  applied  to  the  linear  component  of  the  system  response 
(that  is  to  signals  corresponding  to  such  small  temp  -rature  rises  that  no 
corrections  are  required). 

The  network  connected  between  amplifiers  II  and  III  has  a  transfer 
function  of  the  form  /ITT  This  network  has  the  form  of  alternate  series  resis¬ 
tors  and  shunt  capacitances  (see  Fig.  4)  and  is  a  variant  of  the  standard  so- 
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called  analog  q  network^.  Resistor-capacitor  values  may  be  computed  by  a 
process  involving  continued  polynomial  fractions.  In  practice,  an  algorithm 
due  to  Routh^  is  especially  useful  since  it  is  more  suitable  to  machine 
computation.  A  Fortran  program  for  calculating  the  network  components  for 
a  desired  frequency  range  is  given  in  Appendix  D.  The  measured  response 
of  the  network  of  Fig.  4  approximated  a  '/“s”  function  from  5  Hz  to  40  kHz  to 
better  than  10%.  This  frequency  range  is  adequately  large  to  cover  a  range 
of  test  durations  extending  from  about  2  msec  to  10  msec. 

Amplifier  HI  isolates  the  analog  network  from  loading  effects  of  external 
circuitry  and  provides  a  voltage  gain  to  compensate  for  the  insertion  loss  of 
the  network. 

A  breadboard  model  of  the  complete  compensated  analog  circuit  is  shown 
in  Fig.  5.  A  large  fraction  of  the  bulk  of  the  components  is  seen  to  be  assoc¬ 
iated  with  the  capacitors  of  the  analog  network.  The  operational  amplifiers 
are  of  the  integrated  circuit  type;  very  compact  and  very  inexpensive. 

For  convenience  in  calibrating  the  compensated  analog  network,  parabolic 
type  waveforms  were  applied  as  simulated  signals  from  the  heat  transfer  gage 
to  the  circuit  to  test  its  operation.  With  the  function  generator  disabled,  the 
circuit  produced  a  step  output  for  a  parabolic  input  signal,  as  is  proper.  With 
the  diode  circuit  activated,  and  with  sufficiently  large  amplitude  input  signals, 
the  output  signal,  corresponding  to  heat  flux,  had  the  rising-step  waveform  as 
shown  in  Run  #2  of  Appendix  A.  A  parabolic  signal,  corresponding  to  a  gage- 
measured  temperature  rise  of  1000  UF  in  10  msec  (the  maximum  design 
condition),  was  applied  as  a  real-time  signal  to  the  compensated  analog  net¬ 
work.  The  input  signal  and  the  circuit  output  signal  were  recorded  photo¬ 
graphically  from  an  oscilloscope  display.  The  input  signal  was  manually 
digitized  into  a  temperature -time  tabulation.  These  data  were  then  converted 
to  heat  flux  by  the  digital  computer  program.  Heat  flux  data  scaled  from 
analog  circuit  output  signal  agreed  with  the  computed  data  to  within  a  few 
percent.  Although  time  did  not  permit  extensive  testing  of  the  circuit  using 
a  variety  of  input  waveforms,  the  bandwidth  of  the  circuit  is  sufficiently  broad 
that  all  but  extremely  large  temperature  gradients  should  be  correctly  trans¬ 
lated  into  a  correct  heat  flux. 

Several  observations  on  the  compensated  circuit  are  pertinent  here. 

First  is  the  need  for  the  diode  function  network  to  produce  a  smooth  curve, 
that  is  the  transition  of  conduction  from  one  diode  to  the  next  must  be  carefully 
blended.  Any  perturbations  in  the  function  generator  characteristic  will 
become  painfully  obvious  after  the  differentiating  action  of  the  analog  network 
(see  Appendix  A,  comment  on  Run  #1).  This  problem  was  initially  encountered 
but  it  is  tractable.  Second  is  the  fact  that  two  minor  modifications  would  be 
required  to  the  circuit  of  Fig.  1  to  adapt  it  to  practical  use. 

For  proper  circuit  operation,  the  initial  voltage  drop  across  the  gage 
must  be  nulled  at  the  input  to  amplifier  I  and  the  voltage  input  to  the  diode 
network  must  bear  a  constant  relationship  to  the  gage  temperature  rise.  Since 
the  initial  gage  resistence  and  the  gage  calibration  constant  (change  in  gage 
resistance  for  a  unit  change  in  temperature)  vary  from  gage  to  gage,  calibrated 
controls  can  be  installed  to  vary  gage  current  and  the  gain  of  amplifier  I  to 
make  these  corrections.  Signal-to-noise  ratio  at  low  heating  rates  would  be 
improved  by  a  range  switch.  A  convenient  arrangement  would  use  full-scale 
output  corresponding  to  10,  100,  or  1000  Btu/ft^sec  with  the  diode  function 
generator  inactivated  for  the  lower  two  scales. 
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SECTION  VII 


SUMMARY 


Experimental  studies  were  conducted  to  develop  economically  practical 
methods  of  computing  incident  heat  flux  from  thin-film  resistance  thermometer 
measurements  made  at  hypersonic  shock  tunnel  conditions  that  include  very 
high  heating  rates  and  that  result  in  nonspecific  surface  temperature  histories 
from  the  gage.  Both  digital  and  analog  techniques  that  are  consistent  with 
these  requirements  were  developed.  In  both  methods,  compensation  is  made 
for  thermal  effects  on  the  temperature  coefficient  of  resistance  of  the  thin- 
film  element  as  well  as  on  the  materials  properties  of  the  substrate. 

The  digital  program  accepts  quantized  temperature -time  data  and  pro¬ 
vides  a  highly  accurate  solution  for  the  surface  heat  flux  history  as  a  tabular 
printout.  This  program  was  written  in  Fortran  IV  language  and  made  oper¬ 
ational  on  both  the  IBM  7094  and  IBM  360/65  Model  H  systems.  Since  this 
program  also  evaluates  temperature  and  heat  flux  at  many  planes  within  the 
substrate,  it  is  directly  useful  in  general  studies  of  heat  conduction  in  materials 
of  nonconstant  thermal  properties  at  nonconstant  heat  flux. 

Design,  synthesis,  and  evaluation  of  a  compensated  analog  network  was 
achieved.  This  compact,  solid-state  electronic  circuit  operates  directly,  in 
real  time,  on  the  electrical  output  signal  from  the  thin-film  resistance 
thermometer  to  produce  an  electrical  output  signal  proportional  to  the  true 
heat  flux.  The  circuit  is  simple  in  design,  stable,  and  easy  to  use.  Cost  and 
size  of  the  actual  components  is  such  that  use  of  individual  units  in  each  heat 
transfer  gage  channel  (on  model  tests)  is  economically  feasible. 

Using  the  CAL  analog  computer  facility  (COMCOR),  the  heat  flux  data 
corresponding  to  a  variety  input  temperature -time  functions  was  obtained. 
Generally,  two  solutions  for  each  case  were  obtained,  one  for  constant  sub¬ 
strate  properties  and  one  with  variable  substrate  properties,  to  illustrate 
the  effects  of  these  variables.  In  addition,  these  fundamental  input  functions 
were  also  modified  to  include  tunnel  starting  phenomena,  random  noise,  and 
sinusoidal  modulation.  The  tests  showed  that  measurement  of  the  true  heating 
rate  in  the  steady  flow  portion  of  a  shock  tunnel  test  may  transcend  the  problem 
of  proper  data  reduction.  Cases  were  illustrated  where  the  initial  transient 
portion  of  the  temperature  record  could  seriously  affect  the  heating  rate  even 
at  times  in  excess  of  the  run  duration  of  the  typical  shock  tunnel  because  of 
the  very  long  persistence  (slow  decay)  of  the  transient  effects.  It  is  recom¬ 
mended  that  detailed  studies  be  made  of  tunnel  starting  phenomena  and  initial 
(presteady)  flow  conditions  to  establish  whether  these  factors  do  indeed  con¬ 
found  the  heat  flux  measurement  made  in  the  steady  flow  portions  of  the  test. 
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TABLE  I 


A  Comparison  of  Actual  Temperature  Rise  and  Gage- 
Measured  Temperature  Change  (Relative  to  70  °F) 


Measured  Temperature 
Change,  °F 


Actual  Temperature 
Change*,  °F 


0 


0 


29.9 
127.8 
223.  1 
315.  8 
406.0 
493.5 

578.4 
660.  9 
740.  7 
817.  9 

892.5 


30 

130 

230 

330 

430 

530 

630 

730 

830 

930 

1030 


TABLE  II 


Thermal  Properties  of  #7740  Pyrex  as  a 
Function  of  Temperature 


Temperature 

°F 

c 

Btu  °F" 1  ft'3 

k 

Btu  sec”1  ft" 

70 

25.  03 

2.  17  x  10 

122 

26.  60 

2.  19 

212 

29.  20 

2.  20 

392 

33.  15 

2.  22 

572 

36.  25 

2.  17 

752 

38.  37 

2.  12 

932 

40.  40 

2.06  x  10 

These  values  are  exact  (not  rounded  off). 

>jc* 

These  data  are  taken  from  the  following  source:  Shand,  E.  B.  ,  Glass 
Engineering  Handbook,  McGraw-Hill,  2nd  Ed. ,  1958. 
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TABLE  HI 

Comparison  of  Heat  Flux  Calculations  with  and  without  Corrections 
for  Thermal  Effects  on  Gage  Properties  -- 
Digital  Data  Reduction  Program 


Surface  Temperature*  History  =  7  t}  ^  +  70  °F 


Time,  t, 

Heat  Flux, 

2 

Btu/ft  sec 

sec. 

With  Corrections 

/  =  3.  062  x  103 

7=  4.  593  x  103 

7=  6.  124  x  103 

7=  9.  186  x  103 

.  001 

210.  77 

322. 80 

439.  52 

687. 86 

.  002 

214.42 

331. 35 

455. 24 

721. 03 

.  003 

217. 30 

338.  18 

467.  67 

743.  14 

.  004 

219.  80 

344. 06 

476.  97 

761. 54 

.  005 

222.  07 

349. 28 

484. 57 

777.  81 

.  006 

224.  10 

353. 66 

491. 10 

792.  96 

.  007 

226.01 

357.  32 

497.  06 

807. 28 

.  008 

227.  77 

360. 74 

502.  56 

821. 14 

.  009 

229.48 

363. 78 

507.  83 

834. 54 

.010 

231. 08 

366. 61 

512.  74 

847.  79 

Without  Corrections 

200.00 

300. 00 

400. 00 

600. 00 

The  surface  temperature  history  is  taken  to  be  that  measured  by  the  gage 
(TmofEq.  6). 
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APPENDIX  A 


1.  INTRODUCTION 

Solutions  to  the  nonlinear  heat  diffusion  equation  were  obtained  in  the 
COMCOR  analog  computer  in  terms  of  the  incident  heat  transfer  rates 
corresponding  to  specified  substrate  surface  temperature  histories.  Basically, 
the  primary  input  has  been  a  parabolic  function  with  lesser  emphasis  on  ramp 
and  step  type  inputs.  Perturbations  to  these  basic  functions  were  added  to 
simulate  such  factors  as  noise,  gage  abrasion,  and  tunnel  starting  phenomena. 

Temperature  excursions  of  the  order  of  1000°F  were  used  to  ascertain 
the  full  effects  of  variations  in  substrate  properties  on  the  computed  heat  flux. 
These  effects  are  quantitatively  illustrated  by  obtaining  two  separate  solutions 
for  the  trial  case,  one  with  the  substrate  properties  held  constant  and  one  with 
the  properties  permitted  to  vary  with  temperature.  This  procedure  has  been 
followed  in  all  instances  except  for  those  unique  situations  where  the  principal 
objective  was  to  demonstrate  the  relative  changes  in  heat  flux  resulting  from 
a  modification  in  the  input  function.  In  these  circumstances,  it  was  felt  that 
permitting  the  substrate  properties  to  vary  would  only  add  unnecessary  confusion. 

Smoothing  of  the  data  by  an  electronic  filter  was  provided  as  an  option 
to  demonstrate  the  efficiency  of  this  technique  when  applied  to  various  test 
situations. 

Both  the  input  temperature  function  and  the  resultant  heat  transfer  rate 
were  recorded  by  an  x-y  plotter  on  a  single  graph.  Each  record  is  annotated 
with  a  tabulation  of  the  pertinent  test  variables.  A  total  of  46  trial  solutions 
were  made.  These  are  included  in  this  appendix  and  identified  as  Run  #1 
through  Run  #46.  The  records  are  self-explanatory  to  a  large  degree,  however, 
concise  observations  and  comments  of  the  various  salient  features  of  each 
record  are  given  in  the  commentary  that  follows. 

2.  COMMENTS  ON  ANALOG  COMPUTER  RUNS* 


Run  # 
1 


Test  conditions:  parabolic  temperature  input  with  AT  equal  to 
1000°F  in  10  msec,  initial  temperature  70°F  and  substrate  thermal 
properties  assumed  constant.  Calculated  heating  rate,  q  is 
652.  2  BTU/ft^  sec.  Recorded  trace  shows  slight  "droop"  in  q 
level  with  time  as  a  result  of  the  finite  number  (and  spacing)  of 
substrate  slabs  used  in  programming  the  problem.  Parabolic 
input  is  generated  by  a  diode  function  generator  and  represents  a 
composite  of  individual  straight  line  segments.  Scallop-effect  on 
q  signal  arises  from  the  segmented  nature  of  the  AT  curve. 

Starting  transient  in  q  signal  is  solely  attributable  to  the  recording 
instrument  dynamic  characteristics. 


In  all  cases  the  temperatures  used  correspond  to  actual  substrate  values. 
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Run  # 


2  Test  conditions:  same  as  above  with  the  exception  that  heat 
transfer  gage  substrate  thermal  properties  are  permitted  to  vary 
with  temperature.  Error  in  q  resulting  from  assumption  of 
constant  substrate  properties  is  seen  to  be  negative  and  increases 
in  amplitude  with  increasing  temperature  (at  10  msec  q  is  approxi¬ 
mately  87%  of  the  proper  value).  Most  rapid  change  in  q  (due  to 
changes  in  substrate  properties)  occurs  at  the  early  portion  of 

the  record  (for  AT*Ct^'^)  since  the  rate  of  change  of  temperature 
is  greatest  here. 

3  Test  conditions:  same  as  #1  except  that  the  temperature  input 
signal  has  been  processed  through  a  quadratic  filter  with  a 
corner  frequency  of  5KHz  and  a  damping  factor  of  0.  7.  Effect  of 
the  filter  is  relatively  small  and  is  evidenced  by  a  rounding  off  of 
the  scallop  peaks  in  the  q  trace  and  a  small  time  lag  in  both  traces 
relative  to  both  Run  #1  and  Run  #2. 

4  Test  conditions:  identical  to  above  except  that  gage  substrate 
properties  are  permitted  to  vary  with  temperature.  Comments 
for  Run  #2  and  Run  #3  are  appropriate. 

5  Test  conditions:  same  as  #3  except  for  a  change  of  filter  corner 
frequency  to  2.  5  KHz.  Filter  effects  are  more  pronounced: 
"scallop"  amplitude  well  suppressed,  q  trace  initial  rise  time 
still  slower,  time  lags  in  temperature  and  q  traces  more  obvious. 
Lag  caused  by  filtering  can  be  detected  by  noting  temperature 
amplitude  at  some  fixed  time  (see  tabulation  below): 


time 

temp.  ,  °F 

filter 

0.  5  msec 

297 

out 

0.  5  msec 

287 

5  KHz,  4>  =  0.  7 

0.  5  msec 

275 

2.  5  KHz,  4*0.7 

A  second  q  trace  was  recorded  with  the  filter  damping  factor 
changed  to  1.  0  (critical  damping).  Whereas  the  rise  time**  has 
been  considerably  slowed,  the  time  required  to  settle  at  the 
ultimate  steady- state  amplitude  is  the  same  but  without  the  attendant 
overshoot.  Consequently,  for  all  succeeding  cases  employing 
filtering,  <j>  =  1.0. 


A  quadratic  filter  is  characterized  by  a  transfer  function  of  the  following  form: 

~2 —  where  s  is  the  complex  frequency. 


%/ei  = 


s  +24>w  s+w 
o  o 


ju>,  is  the  corner  frequency,  and  4>  is  the  damping  factor. 

Rise  time,  in  the  sense  used  here,  refers  to  the  time  required  for  the  signal 
to  first  attain  the  amplitude  of  the  final  steady  level. 


34 


Run  # 

6  Test  conditions:  same  as  above  with  substrate  properties  variable 
with  temperature.  Note  increase  in  time  lag  in  temperature  trace 
(relative  to  Run  #5)  resulting  from  the  change  in  the  damping 
factor,  4  ,  from  0.  7  to  1.0  (at  0.  5  msec  temperature  is  now 
approximately  262°F).  Except  for  the  transient  portions  of  the 
record  and  the  presence  of  time  lags,  q  amplitudes  are  not 
affected  by  the  filtering  employed. 

7  Test  conditions:  same  as  Run  #1  with  the  addition  of  random  noise 
to  the  temperature  signal.  Despite  the  random  nature  of  the  noise, 
this  particular  10  msec  segment  of  noise  could  be  faithfully 
reproduced  (in  both  a  time  and  amplitude  relationship).  This 
capability  was  necessitated  since  separate  computer  runs  were 
required  to  acquire  first  the  temperature  trace  and  then  the  q 
trace  on  the  x-y  plotter.  With  the  noise  input  reproducible,  there 
is  a  direct  correspondence  of  the  noise  components  on  the  two 
traces.  The  high  level  of  sensitivity  of  the  q  signal  to  fluctuations 
in  temperature  is  dramatically  illustrated. 

8  Test  condition:  same  as  above  with  the  temperature  signal 
processed  through  a  2.  5  KHz  filter. 

9  Test  conditions:  same  as  above  with  filter  corner  frequency 
changed  to  5  KHz.  The  desirability  of  filtering  the  temperature 
signal  is  amply  demonstrated  by  the  results  shown  in  this  record 
and  the  preceding  one. 

10  Test  conditions:  temperature  input  a  50 °F  step  of  1  msec  duration 
intersecting  a  parabolic  temperature  trace  having  a  vertex  at  70 °F. 
This  input  function  is  designed  to  simulate  a  possible  type  of  tunnel 
starting  phenomena.  With  a  step  temperature  change,  the  corre¬ 
sponding  q  varies  inversely  with  the  square  root  of  time  as  shown 
in  the  first  msec  of  the  record.  The  recovery  of  the  q  trace  after 

1  msec  appears  slow  (despite  the  absence  of  filtering)  because 
the  parabola  is  intersected  at  a  point  (  — '  120°F)  where  its  time 
derivative  is  much  smaller  than  it  is  at  the  vertex  where  it  is 
infinite.  Since  q  is  approximately  related  to  the  temperature 
derivative,  its  rate  of  rise  to  the  final  level  is  less  rapid  than  had 
the  parabolic  function  been  intersected  at  the  vertex.  The  final 
value  of  q  is  not  affected  by  the  presence  of  the  1  msec  duration 
step  (compare  with  Run  #l),  however,  the  time  to  reach  steady 
state  is  approximately  4  msec  compared  to  a  few  tenths  msec  for 
the  equivalent  case  of  Run  #1. 

11  Test  conditions:  same  as  above  with  the  substrate  properties 
permitted  to  vary  with  temperature.  Comments  as  above  apply. 

12  Test  conditions:  same  as  Run  #10  with  the  temperature  signal 
processed  through  a  5  KHz  filter.  The  leading  edge  of  the  step 
input  is  now  rounded  and  the  initial  peak  value  of  q  is  greatly 
reduced.  Time  for  q  to  attain  steady  state  level  is  still  approxi¬ 
mately  4  msec  and  unaffected  by  filtering. 
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Run  # 

13  Test  conditions:  identical  with  above  except  filter  corner  frequency 
changed  to  2.  5  KHz. 

14  Test  conditions:  identical  with  above  with  substrate  properties 
permitted  to  vary  with  temperature. 

15  Test  conditions:  temperature  input  comprises  a  50 °F  step  of 

unlimited  duration  with  a  parabolic  temperature  function  super¬ 
imposed  at  1  msec  (vertex  of  the  parabola  is  at  120°F).  This  test 
is  a  variant  of  the  simulated  tunnel  starting  phenomenon  explored 
in  Run  #10.  The  net  q  in  this  case  is  the  sum  of  heat  fluxes  from 
a  step  temperature  (q  and  a  parabolic  temperature 

variation  (q  =  a  constant).  Since  the  decay  of  q  due  to  the  step  is 
very  slow  with  time,  the  net  q  at  10  msec  is  still  considerably 

in  error*  (672  Btu/  ft^sec  compared  with  652  Btu/ft^sec  without 
the  step). 

16  Test  conditions:  same  as  above  with  the  temperature  signal 
processed  through  a  2.  5  KHz  filter.  Initial  peak  amplitudes  of 
the  transients  are  effectively  reduced,  however,  the  after  portions 
of  the  traces  remain  unaffected  (q  @  10  msec  is  still  ~  672  Btu/ 
ft^sec). 

17  Test  conditions:  identical  with  above  except  for  addition  of  random 
noise. 

18  Test  conditions:  basic  parabolic  temperature  input  as  in  first  test 
(Run  #1  >  except  for  the  addition  of  sinusoidal  noise  (10  KHz). 

19  Test  conditions:  as  above  with  temperature  input  processed 
through  a  2.  5  KHz  filter. 

20  Test  conditions:  identical  with  Run  #18  except  that  substrate 
properties  are  permitted  to  vary  with  temperature. 

21  Test  conditions:  same  as  above  with  temperature  input  processed 
through  a  2.  5  KHz  filter**. 

22  Test  conditions:  temperature  input  consists  of  a  50  °F  temperature 
step  of  unlimited  duration  with  a  parabolic  function  superimposed 
at  time  zero  (vertex  of  the  parabola  is  at  120°F).  This  test  is  a 
variant  of  Run  #15  and  as  in  the  previous  case  q  is  the  sum  of  heat 
fluxes  corresponding  to  the  step  and  parabolic  components  of  the  input. 
Amplitude -time  coordinates  for  Runs  #15  and  #22  are  seen  to  be  very 
similar  except  for  the  initial  portions  of  the  record  (t  £  2  msec). 


It  is  presumed  that  only  the  q  associated  with  the  parabolic  component  of 
the  temperature  function  is  sought. 

Starting  with  this  run  a  new  batch  of  coordinate  paper  was  used  on  the  x-y 
plotter  and  had  somewhat  different  grid  dimensions  from  the  initial  batch. 
Note  that  run  duration  now  appears  to  be  9.  9  msec  rather  than  10  msec. 
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Run  # 


Test  conditions:  identical  with  above  with  input  temperature 
signal  conditioned  by  2.  5  KHz  filter.  Except  for  initial  transient 
amplitude  and  time-  lag,  q  signal  is  the  same  as  for  above  test 


Test  conditions:  identical  with  above  with  the  exception  that  the 
amplitude  of  the  step  component  of  the  input  signal  has  been 
reduced  by  one  half  to  25  °F.  Comparison  of  q  traces  from  Run  23 
and  Run  24  shows  that  the  q  component  resulting  from  the  step 
input  has  been  reduced  by  one  half. 

Test  conditions:  same  as  above  except  for  removal  of  2.  5  KHz 
filter.  As  demonstrated  by  many  previous  test  cases,  the  filter 
is  only  effective  in  modifying  the  early  portions  of  the  tes;  records 
and  is  completely  ineffective  at  later  test  times  (presuming  tran¬ 
sients  occur  only  at  times  near  zero). 

Test  conditions:  parabolic  temperature  input(vertex  at  70  °F) 
with  a  50  °F  step  superimposed  at  2  msec.  This  test  simulates 
the  result  of  a  particle  striking  a  heat  transfer  gage  and  pausing 
a  permanent  change  in  gage  resistance  corresponding  to  a  50  °F 
change  in  surface  temperature.  The  slow  characteristic  decay  in 
q  with  time  (for  step  temperature  input)  produces  a  sizable  error 
in  q  at  9-10  msec. 

Test  conditions:  identical  with  above  except  that  the  input  tempera¬ 
ture  signal  is  conditioned  by  a  2.  5  KHz  filter. 

Test  conditions:  parabolic  temperature  input  (vertex  at  70  °F)  with 
50°F  steps  superimposed  at  2  msec  and  5  msec.  As  in  Run  #26, 
this  test  simulates  gage  damage  by  particle  impact.  The  heat 
flux  components  resulting  from  the  steps  are  additive  so  that  the 
absolute  incremental  error  in  q  is  about  twice  as  large  (at  10  msec) 
as  for  Run  #26. 

Test  conditions:  identical  with  above  except  that  the  temperature 
signal  is  conditioned  by  a  2.  5  KHz  filter. 

Test  conditions:  identical  with  Run  #28  except  that  the  amplitude 
of  the  step  temperature  change  is  25  °F.  Since  q  corresponding 
to  a  step  temperature  changes  extremely  slowly  after  the  initial 
few  milliseconds,  for  the  particular  temperature  record  used  here 
the  net  q  increment  at  10  msec  represents  approximately  equal 
increments  from  each  step.  Thus,  the  q  error  increment  at  10 
msec  (Run  #30)  resulting  from  two  25  °F  steps  is  approximately 
equal  to  that  caused  by  a  single  50  °F  step  (Run  #26). 

Test  conditions:  as  above  with  the  temperature  input  signal 
conditioned  by  a  2.  5  KHz  filter. 
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Run  # 


32 


33 


34 


Test  conditions:  a  50°F  step  input  (with  origin  at  70°F)  conditioned 
by  a  2.  5  KHz  filter.  Addition  of  the  filter  was  necessitated  to 
prevent  the  analog  computer  from  overload  (saturation)  as  a  result 
of  the  step  input  (note  the  change  in  ordinate  scale  for  this  run). 

Test  conditions:  temperature  input  signal  is  a  50  °F  step  (origin 
at  70  °F)  with  a  parabola  (30  °F  rise  in  10  msec*)  superimposed 
starting  at  1  msec.  This  test  is  a  variation  of  Run  #15  with  the 
step  as  the  dominant  component  rather  than  the  parabolic  tempera¬ 
ture  function,  A  condition  is  simulated  wherein  the  starting 
phenomenon  overshadows  the  convective  heating  due  to  the  steady 
flow.  Comparing  Run  #32  and  Run  #33  with  the  net  q  at  10  msec 
is  about  equally  divided  between  the  two  components  (step  and 
parabola).  The  temperature  input  signal  is  conditioned  by  a  2.  5 
KHz  filter. 

Test  conditions:  identical  with  above  except  that  the  parabolic 
component  of  temperature  has  been  changed  to  a  rise  of  60  0  in 
10  msec. 


35  Test  conditions:  parabolic  temperature  input  function  with  a 
horizontal  step  of  0.  5  msec  duration  superimposed  at  2  msec. 

The  characteristic  variation  of  q  with  for  a  constant  tempera¬ 

ture  input  is  seen  in  the  2.  0  to  2.  5  msec  interval.  A  rather  slow 
recovery  of  the  q  trace  follows.  This  slow  recovery  occurs  for 
the  same  reason  given  in  the  comments  for  Run  #10.  The  error 
component  resulting  from  a  horizontal  step  is  seen  to  be  negative 
and  persists  for  the  entire  run  although  decreasing  slowly 
with  time.  The  horizontal  step  is  sometin.es  observed  in  shock 
tunnel  test  records.  Reasons  for  its  occurrence  are  not  known. 

36  Test  conditions:  identical  with  above  except  that  the  temperature 
input  is  conditioned  by  a  2.  5  KHz  filter. 

37  Test  conditions:  parabolic  temperature  input  function  with  two 
horizontal  steps  of  0.  5  msec  duration  superimposed;  one  at  2 
msec  and  the  other  5  msec.  It  is  interesting  to  compare  the  q 
traces  for  Run  #37  with  Run  #29.  In  the  latter  case,  the  pertur¬ 
bation  in  q  corresponding  to  each  of  the  two  vertical  steps  in 
temperature  is  identical  whereas  in  the  former  case  the  q  distur¬ 
bances  corresponding  to  the  two  horizontal  steps  in  the  temperature 
function  are  of  different  magnitude  (although  of  the  same  general 
shape).  The  reasons  for  this  behavior  are  to  be  found  in  the 
comments  for  Run  #10.  As  stated  therein,  to  a  crude  approxi¬ 
mation,  the  q  signal  is  dependent  on  the  slope  (time  derivative) 

of  the  temperature  function.  In  a  region  of  sudden  temperature 
change  (such  as  a  vertical  or  horizontal  step)  the  changes  in  the 


10  msec  as  measured  from  the  origin  of  the  parabola  and  not  from  t  =  0  on 
the  record. 
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Run  # 
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39 

40 

41 


42 


43 


44 


45 

46 


q  signal  will  be  related  to  the  relative  changes  in  the  slopes  of  the 
temperature  signal.  In  the  case  of  Run  #10,  the  slope  of  the 
vertical  step  is  infinite  so  that  the  q  response  is  the  same  regard¬ 
less  of  the  slope  of  the  parabola  function  at  which  the  step  is  applied 
(that  is,  the  q  response  is  dominated  by  the  step  disturbance).  In 
Run  #37,  the  slope  of  the  horizontal  steps  is  zero,  therefore,  the 
slope  of  the  parabolic  function  dominates  the  changes  in  the  q 
signal.  Since  the  slope  of  the  parabola  changes  continuously  with 
time,  it  follows  the  relative  time  position  of  the  horizontal  step 
will  determine  the  corresponding  q  response.  The  relative  effect 
of  a  horizontal  step  on  the  q  signal  decreases  as  the  time  at  which 
the  step  is  applied  increases. 

Test  conditions:  identical  with  above  except  that  tne  temperature 
input  function  is  conditioned  by  a  2.  5  KHz  filter. 

Test  conditions:  parabolic  temperature  input  function  with  a  0.  25 
msec  horizontal  step  at  2  msec. 

Test  conditions:  identical  with  above  except  that  the  temperature 
input  function  is  conditioned  by  a  2.5  KHz  filter. 


Test  conditions:  a  linearly  increasing  temperature  input  function 
with  the  origin  at  70°  and  achieving  an  increase  of  1000°F  in  10 
msj^c.  For  constant  substrate  properties,  the  q  signal  varies  as 

4-  1 


Test  conditions:  identical  with  above  except  substrate  properties 
are  permitted  to  vary  with  temperature.  Because  of  the  linear 
rate  of  increase  of  temperature,  the  q  error  (relative  to  the 
constant  properties  case  above)  increases  continuously. 

Test  conditions:  identical  with  Run  #41  except  that  the  temperature 
function  is  conditioned  by  a  2.5  KHz  filter. 

Test  conditions:  identical  with  above  except  that  substrate  properties 
are  permitted  to  vary  with  temperature. 

Test  conditions:  identical  with  Run  #41  with  random  nqise 
superimposed. 

Test  conditions:  identical  to  above  except  that  the  temperature 
input  signal  is  conditioned  by  a  2.  5  KHz  filter. 
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TEMP- of  HEAT  FLUX  BTU/FT  2  SEC 


TEMP-°F  HEAT  FLUX  BTU/FT2sEC. 
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TEMP-OF  HEAT  FLUX  BTU/ft2seC. 


TEMP- OF  HEAT  FLUX  BTu/FT  2  SEC 


TEMP- °F  HEAT  FLUX  BTU/ft2sEC. 
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EMP-OF  HEAT  FLUX  BTU/ft2sEC. 
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TIME-  MSEC 


TEMP- °F  HEAT  FLUX  BTU/FT 2  SEC 


TIME-  MSEC 
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TEMP- OF  HEAT  FLUX  BTU/FT2SEC. 


I  100 


1000 
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RUN  NO.  18 


INPUT  -  PARABOLA 

CONSTANT  SUBSTRATE  PROPERTIES 

FILTER  -  OUT 

NOISE  -  10  KHz  SINUSOIDAL 


TEMP- OF  HEAT  FLUX  BTU/fT^SEC. 


TIME-  MSEC 
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TEHP-OF  HEAT  FLUX  BTu/ft2seC. 


r 


0)2  34567  89  10 


TIME-  MSEC 
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TEMP- Of  HEAT  FLUX  BTU/ft2sEC. 


TEMP-OF  HEAT  FLUX  BTU/ft2sEC. 


r 


6? 


temp-of  heat  flux  btu/ft  2  SEC 


temp-of  heat  flux  btu/ft  2 SEC. 


TIME-  MSEC 


4-4— -f — |__4_4-  m  no.  29  — l-i-l...;  M4  [  ■ 

INPUT  -  PARABOLA  (VTX  AT  70°)  +  STEPS  (50°  AT  2 MS  AND  5MS)  “t 
CONSTANT  SUBSTRATE  PROPERTIES  itT 

FILTER  -  IN 

FREQUENCY  -  2.5  KHz  - ■ 

DAMP  I  NO  FACTOR  -1.0  T~ 


m 

n 

VI 

n 

m 

TEMP-OF  HEAT  FLUX  BTU/FT  2  SEC 
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TEMP- OF 
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TEMP-OF  HEAT  FLUX  BTU/FT 2 
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TEMP-OF  HEAT  FLUX  BTU/FT2SEC. 


TEHP-OF  HEAT  FLUX  BTU/FT2SEC. 


345  6789  10 


TIME-  MSEC 


TEMP- op  HEAT  FLUX  BTU/FT 2  SEC 


TEMP-OF  HEAT  FLUX  BTU/FT  2  SEC. 
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TEMP-OF  HEAT  FLUX  BTU/FT2sEC. 
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TEMP-Of  HEAT  FLUX  BTU/FT 2  SEC 
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APPENDIX  B 


DIGITAL  COMPUTER  PROGRAM 
W.  D.  Fryer 


The  fundamental  equation  for  one -dimensional  heat  conduction  is 


ST  _  /  <9 

dt  ~  pc(t)  S  X 


With  the  parameters  /o  c(T)  and  k(T)  functions  of  temperature,  the  equation  is 
nonlinear. 


Represent  the  solution  as  T(x,t).  The  objective  is  to  obtain  the  solution 
as  a  function  of  time  given  the  following  initial  conditions 

T  (x,  0)  =  f(x)  =  a  constant,  0  £  x  £  a 

and  boundary  conditions 

T(0,t)  =  T(t),  a  prescribed  surface  temperature  history 
T(a,t)  =  a  constant  end  point. 

We  wish  to  obtain,  T(x,t),  q  Zl  * 

Introduce  the  variable  q. 

ST 

q  h  q(x,t)  =  -k(T)  (Fortran  notation  QD  =  q) 


Substituting  into  the  diffusion  equation,  we  obtain 


<9T  _  f  dj 
dt  ~  SX 


(Fortran  notation  TP  s 


<9T 

ctt 


) 


The  temperature  T  and  the  heat  flow  q  are  evaluated  at  discrete  points 

within  the  substrate  as  illustrated  below.  A  total  of  twenty  temperature  planes 

was  used  (T  ,  T, ,  .  .  .  T  where  n  =  19). 

'  o  1  n 

qos(m,  /) 


Qd(0 


QDM 

— I — 


QD(J) 
— I - 


I 


3 


a. 


»-/ 


■» 


Q. 


Q,  Qs 

QMfiKK(t)  QMAKK(z)  QHARKfa '  QMAKK(H-) 
th 

=  position  of  i  temperature  plane  T(i),  i  =  0,n  (Fortran  symbol  XMARK) 

=  position  of  i-1  heat  flow  plane  QD(i-l)  i  =  0,  n  +  1  (Fortran  symbol  QMARK) 


Since  the  greatest  changes  in  temperature  occur  near  the  surface  ( x —  0), 
the  planes  are  not  spaced  equally.  In  this  program  the  plane  spacing  is  defined 
by  a  cubic  function. 

2  3 

p(z)  =  5z  +  y-  +  ^3- 
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,0' 

1- 


X.  =  XMARK  (i)  =  p(i)  i  =  l.n 

Q.  +  j  -  QMARK  (i  +  1)  =  p(i  -  1/2)  Qj  =  QMARK(l)  =  0 

This  transformation  in  the  spatial  coordinate  provides  a  close  spacing  of  the 
planes  near  the  surface. 


To  calculate  the  partial  derivative  with  respect  to  x,  the  quantity 
is  required. 


<?x/  c?z 


■  4  P('Z)  =  5  +  z  +  z2 

9  z 

At  the  temperature  planes 

f 

DZTT(i)  =  10‘6  (5  +  z  +  z2),  z  =  i,  i  =  1,  n. 

At  the  heat  flow  planes 

DZTQD(i)  =  10'6  (5  +  z  +  z2),  z  =  i  -  1/2,  i  =  l.n. 

A  fp.ctor  of  10  k  is  used  to  obtain  a  proper  scale  factor. 

In  computing  the  transient  solution  it  was  advantageous  to  use  a  change 
in  variable  for  the  time  t.  This  change  was  of  the  form  t  =  T  ^  with  equal 
step."  of  the  variable  X  used  for  integration.  It  was  possible  by  this  method 
to  avoid  the  singularity  near  t  =  0  since  the  time  derivative  is  infinite  at  that 
point  for  an  input  function  of  the  form  T  =  t^'  . 

The  substitution  t  =  X  (1  -  e”  )  was  made  with  /3  a  parameter  that 
may  be  used  to  control  the  crossover  point  (defined  as  the  value  of  /S 
that  makes  /S  Z  =1).  Thus,  =  [  / - r) J  9TB  \9F(t)  dZ 


A  simple  Runge-Kutta  integration  scheme  is  used  to  advance  the  tempera 
ture  history  for  equal  increments  of  the  variable  X  -  The  derivatives  are 
computed  from  first  order  differences  by  the  following  relation. 


q  =  QD(I) 


K*[T(I)  -  T(I-l) 
DZTQD(I) 


I  =  1,  N 


where  K  is  the  value  of  k(T)  evaluated  at  the  temperature 

T  -  T(D  ±  T(I-l) 

2 


by  the  Fortran  routine  KF(T).  This  routine  fits  the  experimental  data  given 
in  Table  II  for  k(T)  with  two  linear  equations. 


k(T)  =  2. 1765  x  10’4  +  0.111111  x  10?  T  (T  5:  400) 
k(T)  =  2.  221  x  10'4  -  0.  2973  x  10‘7  (T  -  400)  (T  >  400). 


For  the  partial  derivative  of  temperature  with  respect  to  time, 

3T  i _  ££  jrt_ 

=  <oc(t)  dx  dr 
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TP(I)  = 


•  Wf QD(I  ♦  1)  -  QD(I) 
[RHOC*DZTT(I)J 


I  =  1,  N- 1 


where  W  is  the  value  of  dt/ dX  evaluated  by  the  Fortran  function  WF(  X  ) 
(WF  =  1  -  e‘4r  +  fiX  ).  And  RHOC  is  the  value  of  ^cfT) 

evaluated  at  the  temperature  T (I)  by  the  Fortran  function  RHOCF(T).  The 
data  of  Table  II  for  p  c(T)  are  fitted  by  two  linear  functions 

RHOCF  =  23.7  +  0.  02452T,  T  £  420 

RHOCF  =  34.0  +  0.01155  (T  -  420),  T  >  420 

The  Runge-Kutta  integration  routine  will  produce  updated  values  of  T(I), 
I  =  1,  N  -  1.  T(N)  is  held  constant  by  the  boundary  condition  while  T(0)  is 
determined  from  the  function  TOF  (  f  )  which  is  used  to  generate  a  parabolic 
surface  temperature  history  as  a  trial  problem.  Setting  KONTRL  =  1,  the 
substrate  properties  are  held  constant  at  their  room  temperature  values 
(k(T)  =  2.  17  x  10-4  and  yOc(T)  =  25.03).  Thus  the  accuracy  of  the  program 
may  be  checked  at  these  specific  conditions  since  the  surface  heat  flux  will  be 
constant  and  of  a  value  calculable  exactly  from  Eq.  (3)*.  Empirical,  or 
arbitrary,  surface  temperature  data  may  be  introduced  by  modifying  this 
routine . 


True  time  is  obtained  from  the  variable  X  by  the  Fortran  function 
TIMEF(  T  )  which  evaluates  the  equation  t  =  Z  (1  -  e~ar  ).  Heat  flow  at 
the  surface  is  obtained  by  linearly  extrapolating  the  values  at  the  first  and 
second  heat  flow  planes.  In  the  program  this  extrapolation  is  accomplished 
in  the  following  manner. 

QDS (M ,  1)  =  QD(  1 )  +  EXTR[QD(2)  -  QD(1)] 


where  EXTR  = 


QMARK(2) 

QMARK(2)  -  QMARK(3) 


In  the  present  status,  the  program  will  accept  the  following  input 
parameters,  if  no  input  parameter  isread  in,  then  the  default  values  are 
automatically  supplied. 


Fortran  Symbol 


Definition 


Default  Value 


BETA 

THICKN 

DTAU 

KONTRL 


TOAMPL 


is  used  in  the  Fortran  function  WF(  Z  ) 
substrate  or  slab  thickness 

ATT  is  an  equally  spaced  time  interval  for 
the  Runge-Kutta  integration 
equal  to  1,  pc  and  k  are  kept  constant; 
greater  than  1,  pc  and  k  vary  with 
temperature 

generates  surface  temperature  input 


surface 


=  TOAMPL*tJ 


5000 

1/16 

12 

5  x  lO"6 


10, 000 


In  the  main  body  of  the  report. 
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C.TITLE:  MAIN  PRPMA#  fop  HEAI-ELQH  PABJJ.AL  JllF-M&EMIIAL -EQUATION 

C  W.  FRVFP  1/67 

COMMON  NM.n'mP.NMM 

COMMON  TAU,0TAU,W,3FTA  _  ..  . .  . . 

COMMON  TO,  T  (10)  ,00(10)  , TP  f  10  )  ,X0  , XMARK  ( 10  )  ,  OM ARK  (30  ) 

COMMON  Q0F(2Olt1O)tTS(20l,30),TIMF(‘>01)  ,OAY(?) 

COMMON  KONTRI. 

COMMON  OZTT(IO)  ,07.100110)  .  ..  .  ._ 

RF  At  KF 

COMMON  /LCTO/  To^MOL 
C0MMON/LCTOP/TPRTFM 

NAMPI  IST/  IND|)T/1FTA,THTCKN,0TAU»K0NTRL,T0AMPLtMMAX 

r 

r  Nm  pfpreffntp  thf  NUMBER  OF  .QO  PLANES  . . 

r 

MM  =  10 
MMP  =  NM+ 1 

NMM  =  NM-1  ... 

XN-=MM 

r 

CALL  DVOCHK 
FALL  OATF(0/\v) 

)N^4  =  o 
C 

r 

r  opp  AULT  VALMrS  FOR  INPUT  OATA 
0  THICKM  is  p  F  AO  IN  FCFT 
pnj  a  =  F  000 . 0 

THTTKN=( 1.0/IA.O) /)?, 0 
OT  AU=.5F-*> 

KONTPl^l  ...  .  . 

TOAMPL  =  1 

MM4X=?oo. 

C  FMO  OF  nFFAMt  T  VA(  IJF  S  FOR  INPUT  OATA 

C 

c 

C  IPRTFMsINTTT  At  PRF-RUN  TjFMPFR  ATIJRF 

R  F  AL  TPPTFM 
I ORTFM=T0 , 0 
r 
r 

r  QKSTFP=NUMHFR  OF  PUNOF-K tJTT A  TNTFGRATTON  STF°S 

TNTFCFR  RXSTFP  _  _ 

pkftfp=?ooo. 

r  .  ...  .......  _ _  _ _ _ _ 

c 

SAVTNT=S4VF  FOR  PRINT  TNTFRVAL  OURING  R-K.  INTEGRATION 
ijo  TO  ?00  VALUCS  OF  TI MF ,QOS, TS, CAN  BE  SAVFO  FOR  PRINT 
TNTFGFR  S  AVT  NT 
SAVINT=10 

0 

c 

1  RF  AO  (S, INPUT)  _ 

C=(THTCKN/(  F.*XN  +  XN**?/?.(HXN**1/1.0)  )*i.0F«-6 

WR I TF ( 6 , 1  ) 

1  FORMAT! 1 H) ) 

WRITF  f6, INPUT) 
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?e(KnNm.r,T.n)  r,n  rn  4 

stop 


4  continuf 
vo  *  o.o 
no  ?  *=],vm 

7  =  T 

7  =  7-0. 5 

1  QMAOIMfnt=(,i.O*Z*0,q*7**?4.7**VA.P»*r. 

owark  m  -  r.r » 

1cfTNF4.0T.oi  nn  TO  POl 

or?  i  NT  Rnn,v^fyo,(ymi>Kin  ,OMAP*(  !♦!)  ,  T»l  ,NM| 

POO  format  ( l  HI  OY  -»PHP|  ANF  infATfONS  IN  MTCPOFFF  T  /// 

A  10X  ltHTFMOCOATU«F  nY  4WQnnT  //  MY  Flfl.7,  F17.?M 
TN04  =  1 
f.O  Tp  qo? 
qni  rr"\|TfMiJF 
pnp  f ONTTNUF 
r 

T O- j  pp  TP  M 

no  10  |al  .  aim 
1  o  T  f  f  )  =f ORTCM 
OP  IP  T  =  1  ,  IVJM 
7  =  I 

OTTT(T  )=f*1  .or-p*(F .0+7 +7**7) 

7  -  7-0. P 

1  »  nr  TOP  (  n  =  f*1  .  pp-p*(  p.o+7+7  **7  1 
TAM  -  0.0 

CALL  9  NOKT I  t  ,MM«,  T  AIJ,OTAU,Tt  TP) 

FYTR  a  QMA°K(n/  (QMAPK  CM-0MA9K  CO)  1 
.1  =  o 

m  -  n 

no  1  on  N= 1 .OYPTPo 
J  =  .1  +  1 

Ml  l  ONOYTf  ?,  NMM,  TAII.OTAIJ,  T,TD) 
pp  Tr(  .I.MF.S*\MMT)  OO  TO  100 
.1  =  0 
M  -  M+1 

no  on  NNs'^.N'40 
Oop ( M, NN 1  =  pn (NN-1 1 
on  T p (  M,  l\l\|  1  a  T  (  NN-  1  5 

nc;  pnp(Mfn  =  AIM!  +  FYT9*fOP(  71-00(1  n 
TPf  M, 1  )  =  to 

TlwrfM)  =  TI mcf ( TAH)*1 OOP. o 

r 

ion  rONTTNlJF 
r 

r  o  A|  nil  ATTONS  P0M°1  FTP  , 

r  opmpdatp  OUTPUT 

po  MNT  1  M')c 

oojmt  l?0,  nAV 

1  ?o  F0PMAT1  1  H  1  ,11V,  qpwpijRrACF  t^MPPPATUPF  ANO  COOT  H  fPTORV,  1 [OY,  7A4/  ) 

UP  1 TF (  p ,  1  00) 

170  TnpMATnY,4HTTMF,QYt1HT,11  Y,4H0nnT/->Y,f,H  (MPFC)  1 
PM  NT  1  40,  (  TTMF(  ()  ♦  TP(T  ,11  » OOP  (  T  ,  1  )  ,1  =  1  ,  MM  AY  ) 

140  F0RMAT(1Y,FP.7,F11 ,7,F11 .?) 

Ipp  r.o  to  i 

TNO 
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C TITLE :  SUBROUTINE  nAUX  FOR  HEAT-FLOW  EQUATION.  .W.QF  3/67 
SURROUT IMP  nMiy 
COMMON  NN.NMP.NMM 
COMMON  TAM,OTAM,W,MFTA 

COMMON  TO,T  (  30)  ,00(30)  ,  TPMO)  ,XOf  XMARK(30)  ,QMARK ( 30) 
COMMON  OOF  M01 , 30)  ,TS(  ?0)  ,  30  »,  TTMF(  301)  ,OAV(  3) 

COMMON  KOMTOf. 

COMMON  07TT (30) f nZTOni  30) 

RFAL  KfKr 

TO  •=  TOP(TMl) 

TMyn  =  ( T(i ) +T0  )  /7.o 
K  =  KF( TMTO) 

00(1)  =  -**(  TH  )-TO)/DZTQO(l) 
no  ! o  (■=•>, MM 
TMTO  ={T(T)  +  T(T-l)  )/■».  o 
K  =  KP(TMTO) 

n  QO ( I )  =  -K*(T( T ) — T ( I  —  1 ))/nZTOO(  I ) 

r 

no  ?o  t  •=  i ,  MMM 
R MOO  =  RHOC c ( T ( T  )  ) 

W  =  VF(TA(J) 

70  TP(I)  =  -M*( on( T+1 )-gn(! ) ) /  (R40C*n7TT( T) ) 
r 

DPTIIRM 

PNO 


CT  (T|  p*  Fliwr  T  T  ON  TOP  (  T  All)  -  PARAMO!  A  WOF  3/ft7 

r 

FMNr T TON  TOP(TAII) 
opa!  (PPTFm 

r  OMMON  NM,NMD,NMM 
TOMMON  n||M  t  OTAIJ  ,W,3PTA 

r OMMON  TO, T( 30) ,00(30) f TP( 30 ) ,X0, KMAPK (30) f OMARK ( 30) 
r OMMON  ons ( 70) , 30) ,TS( 701 , 30 ) ,T( MF (70) ) ,OAV( 7) 

common  /|  r to /  toamrl 
COMMON/!  rTOP/TPRTFM 
TT  M  =  TI  MPP  (  T  AM) 
y  =  TOAMPl  *SO°T  (  TTM) 

TOPA  =  (1.0— S0RT(1. 0-5. )RF-4*T))/3. MRP-4 
TOp  =  Toc A  +  T  ootpm 
OF TURN 
PMO 
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r.  T  T  Tl  «= .  H)Mf  T  I OM  *c(T).  CONSTANT 
OPAL  P'INrTfnN  KFfT  I 

COMMON  !0»|0°C  (1  24??  ) 

COMMON  K°MTf?| 

!c ( KONTO|  . OT . I )  00  TO  »n 
*c  - 
RPT(IOK) 

r 

00  I c( T.OT ,400, o )  on  fo  70 

KC  =  ?.)7a<;c_4  f  0.niH1F-7*T 
OPTUPW 

™  «r.  =  ,  _  o. ->07?p-7*(  r-40?  .n  ) 

periJRv 

TNO 


WOF 


rTTT|F:  RHnrc  frn^STAMT)  v, 

P'JNJrTjoKj  DMnrrfr) 
fOM M 0 M  TOMqor  n  ->47->  J 

rnMMPN  K  0*'Tf>| 

T^ft'OMrqi  .OT.I  )  on  rn  ->o 
PMnrc  -  ">s.n7 
O^TIJRM 
r 

90  I  c  (  T,  OT  ,4?0,  n  j  on  TO  ?o 
OMnrn-  -  77,7  n,  07457*7 

RPTIIPM 

to  P(-iorr  =  v,,o  +  o.ni  !*;*i*f  T-4?n,oj 
RCTKOM 
PMO 


7/47 
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CTTTl.F:  ruNCTinN  TIMPF(TAU)  WO®  1/67 

FUNCTTON  TTMFF(TA(M 
COMMON  NM» NMO» NMM 
COMMON  nuM*OTAU,W,RFTA 

COMMON  TO,T(in),0O(in)  ,TP(10)  ,X0,X*ARK(70)  .OMARKCAOI 
TOMMON  OOF (?01,7O),TS(?0!,10),TtMF(701  ),nAV(?1 
in  ARO  =  RFTAotaij 

fP(  Apr,. |  T.7.0F-7)  qn  TO  10 

I *=  f  apo.ot.io.o)  r,n  rq  40 
90  ttmff  =  tau*( 1 .o-pyp(-aroi) 

PFTtJPM 

in  Tf  Mrp  *  TAU*AO0 
pc  TURN 

t,n  TTMFF  =  T  AH 
PCTu°N 
tno 


CTTTIFj  PIINOTTOM  WFITAtn  V'OF  1/47 

r 

FIINCT  I  ON  WC(TAU) 

common  mm,mmp,nmm 
COMMON  ni|M,  HTA'I,  W,  RFTA 

rOMMDN  T0,T(1O)  ,00(10)  ,  TP  (10)  ,XO,YM  ARK  ( 10 ) , OMAR K ( 10 ) 
'■OMMOM  qn<;(?ni,  loi  ,TS(?oi,10),TTMF(?011,OAV(?1 

r. 

AR  0  =  RFTA+TAIJ 
TF(  AR0.0T.7R.nj  co  TO  ’0 
JC(  A»0.|  T.  1  .OP-7)  00  TO  10 
F  =  FX  °  (  -  ARn  ) 

(JP  -  1.0-F  P*ARO 
RCTMPN 
in  uf  ■=  l  ,  n 
PPTURN 

in  =  ?.n*AQO 
PPTUPN 
FNO 
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.  * 1 1)  w  r 


1 1 


•  **  •  WW& 09!  v^WPt^ 


(IBM  t*o> 


rfncT  r  rnHKT  STNOie  PBFFTS'^N  BIJNCr  KUTT  ^  (TRM 

r 

r  ***  *********  *****  ***  *  **  ****  ******  ************************************** 


SURPnilTINF  RNOKT 
r M  on*-> 

niiponSF 

SIM'', |F  pppMSTON  PIJNSF-KUTTA  INTFOPATTON  (HUM'S  MOT!  Pff  M- 
T  HMI  tno  m  (M=1  TO  TO)  nTFFcPFNTI  M.  FOUAMnNS.  FAFH  TAIL 
nr  OM0KT  f  A I  ISP  S  fMTPO RATION  CH0  f)NF  STFO. 

MOTF  ,  All  |  B|  F  fOFrTSWM  IS  ')Scn  IN  I MF  P  F  MFNT  I  NH  TWOc- 
PFNnFMT  VAR  I  AP|  F. 


•  is  a  cr 

r  ai  i 


RNSKT  ( TNI T ,M, Y,MtY,vo} 


ncSFRlOTION  OP  PAPAMFTFRS 

T  AIT  T  =1  POP  INTTJ  AI.T7ATION, 

r->  p OR  TMTFr,f>  A  Te  ONF  STFO 

m  mijmtfp  of  nosT-nonPR  0  TPFPppNT  I  AI.  FQiJATTONS 
y  INnroPNnFMT  vapiaimf  (STNOI.f  PRFCTSfnN) 

M  STPR  S  T  7  F  TSIMOI.P  PPFFTSTnN) 

v  ARRAY  OF  OFPFNTFMT  VARTARIFS  (ST NOT  F  POPflSinN) 

vn  A  P  p  A  V  n  P  OFP  I  V AT IVF  S  OP  Y  MPT  Y  (SINP.IF  PPFf  I  ST  ON) 


PP YAPKS 
i  ISP  P 


ilSPP  MOST  SUPPLY  A  S'lRonilTTNP  OAUY,  Sljrw  THAT  'FAIL  OAKY* 
(VITHO'IT  FYPIIPTT  AROMMPNTS)  CA'ISFS  T HP  YO  VAMJPS  TO  PF 
r A|  PIILATFO,  IJSTNP,  Y  ANn  VALIJFS  FROM  thf  y  AP°AY. 


TNTTi  Al  T  7ATjnN  (TMJT  =  1)  TS  PFO'ITPFO  AFT«  WSFP  SFTS 
t»'TT|A|  VA|'|PS  OF  Y  ANO  TWF  y-APPAY.  p  c_  f  Mf  t  T  AL  !  7  AT  f  ON 
TP  MOT  MPPPSSAPV  pno  f HANrr S  TN  STFP  ST’P,  M. 


P'JOpni|TTNPS  ANT  PIJNPTTHN  Si  IB  PROS  PA  MS 
O  A 1 1 V  (SCC  PPFVin'IS  RPMAPYS) 


PFO'ITPFO 


P*  MeTMnn  * 

p*  THP  F.Y.  B|||M  «nn  I F  I P.  AT  T  ON  Oc  TWP  PIINOF-YIITTA  PHIJPTH-  * 

C*  opofr  wrrMnn  ts  IISFO.  SFF  * 

P*  iMATHrMATTPS  OF  r  r|MO||T  ATT  HN '  »  APPTI,  PA?  PP  17/S-lqT  * 

f*  * 

r-H  A'ITmoo  n a t p  * 

r  a  fp yfp  * 

C*  PHPNFLL  ,CBn,  I  AB  * 

F*  OPTOPPP,  (OlSiA  * 

f*  * 

f*  nptf,  rnni^'F,  is  AnAOTFn  from  somapf  partff  * 

F*  TMTCCPATOP  OROORAMS  PTNT1  ANO  OfNT.  * 

r*  OFfYS  ripnso,  HIOnSA.  * 

r+  * 

T***-****'**************************************************************** 

r 

SiRPOUTTMC  PMSYT  (  T  NT  T  v  Nf  /Y  /  ,  H,  V  t  Y"  ) 


0  P  Al  Y  »  H,  Y  (  T  )»YP(I) 

p  p  Al  P  ( in  | ,  n(oni 
nni|P|  p  ppfCTSTON  YF>f  ho  94 


TNTT 


C 

r,n  to  (4,  PI , 

4  rau  naux 

5  Y0  =  Y 

PPTUPN 

r 

0  HP  =  H 

on  in  f=l,M 
PIT I  *  H*vo( j ) 

vm  =  YHI  ♦  o.  •>*»!!» 
in  cm  =  PI1 1 

c 

x  =  xn  ♦  o.*>no*Hn 
ran  naux 
c 

nn  ■Jn  !*l,M 

PHI  =  H*VPU> 

vm  =  veil  ♦  a.wiPiiHim 
Q(  H  =  Q(  H/4.P 
r 

TALL  PAIJX  . 

r 

nn  on  t  *  i,N 

PIT)  -  H*VPH  )  -  0 ,  <5*P  (  T  ) 

vm  =  vn  i  ♦  p ( t  i 
pp  QH)  =  om  -  ofT) 
r 

xn  =  vn  «.  hp 
x  -  xn 
rail  naijv 
r 

nn  40  T  =  1 »N 

P(  T 1  =  H*vpm  f  (°UHP(TH 
40  YU)  =  V  f  T  I  ♦  fnm+P(H/4.0) 
r 


rau.  naux 

rptiipn 

FNO 
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APPENDIX  C 


DERIVATION  OF  CORRECTION  FUNCTIONS  FOR  THE 
COMPENSATED  ANALOG  CIRCUIT 

S.  A.  R.  Ayyar 


This  appendix  is  comprised  of  three  sections.  Part  I  is  a  summary 
while  Part  II  discusses  the  mathematics  associated  with  the  correction  for 
the  nonlinear  resistance -temperature  function  of  the  gage  resistance  element 
and  Part  III  presents  the  mathematics  underlying  the  compensation  for  the 
variation  of  substrate  thermal  properties  with  temperature. 

I.  SUMMARY 

As  the  thin-film  resistance  thermometer  is  used  to  measure  high  heat 
flux,  the  surface  temperature  may  be  carried  into  the  range  where  the  thermal 
coefficient  of  resistivity  of  the  resistance  element  and  the  thermal  conductivity 
and  diffusivity  of  the  substrate  material  are  affected.  Recent  research  has 
shown  that  these  temperature  sensitive  effects  can  be  combined  mathematically 
into  a  single  correction  factor  and  applied  as  an  operation  on  the  gage  electri¬ 
cal  output.  This  means,  in  practice,  the  gage  output  voltage  is  subjected  to 
electrical  compensation  before  it  is  introduced  as  the  input  to  a  conventional 
heat  transfer  analog  network  to  provide  a  signal  proportional  to  the  true  heat 
flux. 


The  correction  factor,  which  is  a  curve  with  voltage  output  as  abscissa 
and  transformed  temperature  as  the  ordinate,  is  the  result  of  two  effects. 
First,  we  establish  the  relationship  between  voltage  output  and  temperature 
rise  of  the  gage  considering  the  thermal  coefficient  of  resistivity  of  the  resis¬ 
tance  element.  This  relationship  for  a  particular  gage  under  study  takes  the 
following  form: 

=  0.  386  x  104  -  V  0.  149  x  108  -  0.  0771  x  108  (C-)) 

§o 

QJF  =  temperature  rise,  °F 
uq  =  gage  voltage  rise 

U  =  initial  gage  voltage 
60 

The  derivation  of  this  expression  is  shown  in  Part  II. 

Next,  we  establish  the  relationship  between  temperature  rise  obtained 
above  and  transformed  temperature  rise  which  takes  into  consideration  the 
thermal  conductivity  and  diffusivity  effects  of  the  substrate  material.  This 
relationship  takes  the  following  form: 

*£  =  e»[*-(n°-K ]  <c-2) 


=  transformed  temperature  rise  in  °C 
=  temperature  rise  in  °C 
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For  the  particular  gage  and  substrate  thermal  properties  provided,  the 
following  value  is  obtained  for  the  constant  a: 

a  =  -,  00096/  °C  rise 

Since  we  have  6^  in  °F  from  Eq.  (C-l)*,  the  relationship  (C-2)  becomes 

Jo)  - 

*ux  =  (°F)  [1  +  .0001882  (|)  (eF)]  (C-3) 

Here,  (o) 

0 0  =  transformed  temperature  rise  in  °F 
6 -  temperature  rise  in  °F. 

The  derivation  of  Eq.  (C-3)  is  shown  in  Part  III. 

II.  EFFECTS  OF  GAGE  NONLINEAR  RESISTANCE  VARIATION 

A.  Mathematics 

The  resistance-temperature  function  of  the  thin-film  resistance 
element  is  defined  by  the  second  degree  equation.  ^ 

R™  _  ~  =  a  +b  T-c  T2 
T  -  T  o  o  o 

=  a  +  a2  (T  -  70°)  -  &3  (T  -  70°)2  (C-4) 

=  value  of  gage  resistance  at  T  =  70  °F. 

Differentiating  Eq.  (C-4) 
dR 

=  0  +  a2  -  2  a3  (T  -  70°)  (C-5) 

The  constant  a^  is  the  value  of  gage  constant  k  at  T  =  70  °F,  i.  e.  , 
dR  _ 

"  dTT  =  70°F  "  T  =  70° 

Substituting  into  Eq.  (C-5)  and  using  experimentally  determined  values  for 
the  coefficients  ,  Eq.  (C-6)  is  obtained. 


dJl  .  4L  (r-ro-)] 

dTT-T  dTT=10*F  L  *T;70 


(C-6) 


r 

rr-  70  ‘F 


/  -  2.5  7  *  15  ( T-70\ 


The  temperature  coefficient  of  resistance  of  the  gage,  ^  ,  may  be  defined  as 
follows . 


1 


’  Eq.  (C-2)  was  derived  in  the  Centigrade  unit  to  afford  a  direct  comparison 
with  information  given  in  the  referenced  material. 
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*9a  <*e 


/>T"-T  ''  *»„ 


^Tm  70  V  "  5~ 


*6  \  dQ  /  r,7(3V 


In  Ref.  7  the  temperature  coefficient  of  resistance  is  shown  to  vary 
according  to  the  following  function 


fi(*)  =  /*,('-  *o) 


By  integration 


a 

f  fa-to?6)*6  ^  (Our-J  9tr) 


(C-7) 


Thus 


R  q  “  Af\a  A  A?  /  2.  \ 

-hr3'-  ir  =  e°(e"  -  r 


(C-8) 


Rearranging  the  terms  of  Eq.  (C-8),  we  obtain 

AR 


2  ‘  6<*  + 

3o 


=  0 


Solving  for  0 


9 


ur 


=  x  -  /XT 

*  V 


^AR 

0O  *RS 


(C-9) 


Using  the  data  given  in  Ref.  7,  ?  =  2.  59  x  10'^/  °F  and  ^  =  1  x  10”^/ 

JF,  and  substituting  into  Eq.  (C-9) 


0~  = 


/ _ 

2.57*  /O’*- 


/— L-  - 

y  (z. si x to-*)3- 


f 


2AR _ 


(C- 10) 
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Notation 


R 

T 


e 


^ur 

7 

A  R  = 


u 


o 


res  istance 

conditions  at  temperature  T 

gage  resistance  at  0°C 

initial  resistance  of  the  thin  film 

temperature  coefficient  of  resistance  of  the  thin  film 

temperature  rise  measured  from  its  initial  level 

gage  suriace  temperature  rise 

constant  as  defined  in  the  equation  (6) 

change  in  resistance  of  the  thin  film 

voltage  rise  in  the  gage 

initial  voltage  across  the  gage. 


B.  Practical  Data 


By  maintaining  a  constant  excitation  current  through  the  resistance 
element  of  the  gage,  the  following  relation  will  hold 

A  /? 


Substituting  this  relation  in  Eq.  (C-iO) 


eurM  - 


I 


2  U, 


2.5<fx  IO -* 


(2.59  *,<,-+)* 


(  d 6  1 

3 O  '  'T=  7o*F 


?x  td*)U3  (C-ll) 


Equation  (C- 11)  is  the  general  relationship.  For  circuit  design  purposes , 
some  typical  data  were  substituted  for  gage  parameters 

ff  =  loo  ohms  f  =  o.i  ohn/‘F 

yo  \  dt  /Tsyo.F  / 

_ _ - _ _  (C-12) 

&u r  (°F )  ~  0.38  A.  x  10*  —  y  O.  Id-?  x  IO8  -  O  .07  7  /  *  /  O*  — 

U3o 

HI.  EFFECTS  OF  SUBSTRATE  PROPERTY  VARIATION 
A.  Mathematics 


The  derivation  of  the  correction  formula,  Eq.  (C-2),  is  taken 
from  Ref.  14  except  for  some  changes  of  notation.  The  equation  governing 
the  thermal  diffusion  is 


/oc  (e) 


do 


6 


(C-13) 


Here , 

P  -  density 

c  =  specific  heat 

&  -  temperature  rise  above  initial  temperature 

t  =  time 

y  =  distance  below  surface  of  slab 

k  =  thermal  conductivity 


(C-  14) 


We  introduce  (f>  as  the  transformed  temperature 

9 


/ 


0  =  r 


Here  kQ  is  the  thermal  conductivity  at  the  initial  temperature.  In  Ref.  14,  a 
constant  b  is  introduced  which  reflects  a  linear  relationship  between  k  and  kQ. 
In  our  case,  the  k  variation  with  temperature  in  the  range  under  consideration 
is  such  as  to  make  the  relationship  0  versus  9  almost  equal  (see  Table  II). 


to  Eq.  (13). 


Proceeding  as  per  Ref.  14,  transformation  of  Eq.  (14)  is  applied 


S<P 

7t  -  *  77*  =  ° 


(C-  15) 


which  is  the  form  in  which  diffusion  equation  is  customarily  given.  Equation 
(C-15)  is  nonlinear  since  the  diffusivity  K  is  a  function  of  temperature,  i.  e.  , 

K  =  K  (  (0  )  (C- 16) 

Hartunian  and  Varwig^  applied  a  perturbation  technique  to  linearize  Eq.  (C-15). 


Here  we  take 


K  =  K(0)  +  K*1 ) 


(C-17) 


Equation  (C-17)  is  substituted  into  Eq.  (C-15)  and  terms  of  like  order  collected 
to  provide  the  zeroth  order  equation 

d  <t>‘}  (o)  S  *  (a) 

-  K  -.—ST-  -  O  (C- 18) 


St 


Sc* 


and  the  first  order  equation 


d  <t)0>  _  ,  (.) 

St  c?y  * 


K 


0) 


dtp 


(°) 


K  (°)  St 


(C-  19) 


by 


The  heat  flux  q  at  the  boundary  of  the  semi-infinite  body  is  given 

A  -  _  jl  S(t>  I 
*  -  *•  (C-20) 

where  the  subscript  w  refers  to  the  wall  surface,  y  =  0. 


The  solution  to  the  zero^  order  equation  is 

f-% 


4,<‘>  -  2  h 


Kla)t 


(C-21) 


The  first  order  solution  obtained  by  Hartunian  and  Varwig  is  based 
on  taking  k(^)  equal  to  the  constant  Kq,  which  is  the  diffusivity  at  initial 
conditions  and  then  further  taking  proportional  to  i.e., 

(o)  (C-22) 


K 0)  =  cl  Kg  <t> 

With  this  relation,  the  first  order  solution  then  provides 

0)  /  a  * 


0, 


W 


% 


(C-23) 
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Using  Eq,  (C-21 ) 


with  which 


(C-24) 


K 


(C-25) 


/  \  ^  (o) 

The  particular  solution  required  is  that  for  0^  .  Assuming  yy  a-^>CtT<<  I  > 

Eq.  (C-25)  may  be  inverted  without  difficulty  to  yield  a  polynomial  in  . 


( C  -  2  6) 


Using  the  data  on  substrate  properties  given  in  Table  II,  the 
diffusivity  was  approximated  by  the  linear  function  K  =  KQ  (1  +  a  (p  )  v/ith  "a" 
evaluated  from  the  slope  of  K  versus  <p  .  Substituting  the  value  a  =  -0.00096/ 
°C  into  Eq.  (C-26)  and  converting  into  Fahrenheit,  the  desired  correction 
function  is  obtained 

0.  (°F  .  )=  1(°F  .  )  [1  +.  0001882  (5/9)  0/jr(°F  .  )]  (C-27) 

(»  rise  to'  rise  L  tor'  rise  J 
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APPENDIX  D 


COMPUTER  PROGRAM  FOR  DEFINING  COMPONENT  VALUES 
FOR  THE  ANALOG  NETWORK 


The  program  for  designing  the  analog  network  having  a  frequency  response 
k  VT  is  written  in  Fortran  and  was  run  on  a  Quicktran  terminal  served  by  an 
IBM  7044  computer.  With  minor  modifications,  it  should  run  on  any  IBM 
computer  as  well  as  several  other  machines.  Input  variables  are  written  into 
the  program  since  it  is  not  anticipated  that  several  networks  would  be  designed 
at  any  one  time.  These  variables  specify  the  lowest  frequency  at  which  the 
response  of  the  network  approximates  (see  Note  1)  and  the  range  of 

frequency  coverage,  specified  as  the  logarithm  of  the  ratio  of  the  highest 
usable  frequency  to  the  lowest  usable  frequency  (see  Note  2).  There  are 
eight  RC  sections  in  the  network.  According  to  Skinner^,  a  network  of  n  sec¬ 
tions  should  adequately  cover  (2n  -  1)  octaves,  i.  e.  ,  an  eight  section  network 
is  useful  over  15  octaves  or  4.  5  decades. 

Output  is  specified  in  the  form  of  resistance  (ohmsl  and  capacitance 
(farads)  values  of  the  elements  of  the  network.  In  the  program  printout,  the 
first  value  printed  is  a  zero’1'  succeeded  by  a  one  ohm  terminating  resistor, 
followed  by  a  shunt  capacitance,  a  series  resistor,  and  so  on  for  a  total  of 
seventeen  elements.  The  values  are  scaled  to  suit  the  designer  by  multiplying 
resistances  by  a  convenient  factor  and  dividing  capacitances  by  the  same 
factor.  A  listing  of  the  program  follows: 

LIST 

PROGRAM  POLY 
DIMENSION  DUP(1  70),  AM(18) 

DIMENSION  AR(9,  18),  ELEMENT(18) 

DIMENSION  A(18),  D(18),  P(*),  Z(*) 

INTEGER  AO,  DO,  PO,  ZO 
WRITE  (6,  92) 

92  FORMAT  (8X,  8HSTARTING) 

AO  =  0 
PO  =  1 
DO  =  0 
ZO  =  1 
NP  =  1 
NZ  =  1 

DO  10  J  =  l,  18 
ELEMNT(J)  =  0. 

DO  10  1=1,  9 
A(J)=0. 

D(J)  =  0. 

10  AR(I,  J)=0. 

A ( 1 )=  1 . 

D(  1 )-  1 .  ...(1) 


This  first  value  has  no  significance  and  should  be  disregarded. 
*  Lower  limit  (Hz). 
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Pl=. 003 

DO  131  J  -  1 ,  A  0 

P1=P1*6.  28318 

...(2) 

I=AO- J  + 1 

AN  4.  5 

131  AR(I,  1  )  +  D(J) 

M  =  8 

DO  132  J=l,  DO 

M2  =  2*M 

I=DO- J  + 1 

RD  M2  + 1 

132  AR(I,  2)=A(J) 

DO  1  1=1 ,  M2 

91  FORMAT(8X,  218) 

AI  I 

130  CONTINUE 

RN  AI*AN 

DO  20  J  =  3,  18 

R-RN/RD 

DO  20  1=1, 8 

1  AM(I)  =  10.  *#R 

IF(AR(1 ,  J-  1 ) ) 2 5 ,  25,  19 

DO  2  1=1,  M 

19  CONTINUE 

L=2*I=1 

AR(I,  J)=AR(I+1,  J-2)-AR(l,  J-2)* 

Z(I)-P1*AM(L) 

AR(I  +  1,  J-1)/AR(1,  J-l) 

2  P(I)  =  P1*AM(L+1) 

20  CONTINUE 

WRITE  (6,  99)P 

...(3) 

25  CONTINUE 

WRITE  (6,  99)Z 

•  •  • (4) 

DO  30  L=2,  18 

NP=M 

IF(AR(  1 ,  L))3 1 ,  3 1 ,  29 

NZ=M 

29  ELEMNT(L)=AR(  1 ,  L  -  1 )/  AR(  1 ,  L) 

AO-1 

30  CONTINUE 

DO=  1 

31  CONTINUE 

1=  1 

WRITE  (6,  90)AR  ...(7) 

DO  111  N=l,  NZ 

90  FORMAT(lX,  9E12.  4) 

DO  110  IND=1,  AO 

WRITE  (6,  90)ELEMNT  ...  (8) 

I=AO-IND+l 

STOP 

A(I+1)=A(I+1)*Z(N)+A(I) 

END 

110  CONTINUE 

A(1)=A(1)*Z(N) 

AO-AO+1 

111  CONTINUE 

D0121  N=l,  NP 

DO120  IND=1,  DO 

I=DO-IND+l 

D(I+1)  =  D(I  +  1)*P(N)+D(I) 

120  CONTINUE 

D(1)=D(1)*P(N) 

DO-DO+1 

121  CONTINUE 

WRITE  (6,  99)A 

WRITE  (6,  99)D 

...(6) 

99  FORMAT(lX,  9E12.  4) 

^Logarithmic  range 

3Poles 

47 

Zeros 

Polynomial  (generated  from 

zeros ) 

Polynomial  (generated  from  poles) 

7 

Array  for  calculating  with  Routh's 

algorithm  (Ref.  15) 

Element  values 
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